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Folding of RNA is subject to a competition between entropy, relevant at high 
temperatures, and the random, or random looking, sequence, determining the low- 
temperature phase. It is known from numerical simulations that for random as 
well as biological sequences, high- and low-temperature phases are different, e.g. the 
exponent p describing the pairing probability between two bases is p = | in the 
high-temperature phase, and p ~ | in the low-temperature (glass) phase. Here, we 
present, for random sequences, a field theory of the phase transition separating high- 
and low-temperature phases. We establish the existence of the latter by showing that 
the underlying theory is renormalizable to all orders in perturbation theory. We test 
this result via an explicit 2-loop calculation, which yields p « 1.36 at the transition, 
as well as diverse other critical exponents, including the response to an applied 
external force (denaturation transition). 



I. INTRODUCTION 
A. Random RNA 

Together with DNA and proteins, RNA plays a key role in biology. As such it is important 
to understand its spatial conformations. While for protein the lowest-energy fold depends 
strongly on the chemical constitution, and is only tractable numerically, the problem for 
RNA is simpler, due to a clear separation in energy scales between primary structure (the 
sequence), secondary structure (pairing of bases in a fold) and tertiary structure (embedding 
of a fold in 3-d space). 

RNA molecules consist of 4 bases - adenine, guanine, cytidine and uracil - which are 
attached to a flexible sugar-phosphate backbone. In contrast to duplex DNA molecules 
(where uracil is replaced with thymine), there does not exist an independent complementary 
strand, and the RNA molecule folds back on itself. Experimentally important (see e.g. pQ) is 
the observation, that topologically intertwined pairings such as knots and pseudo-knots do 
not seem to play a crucial role for the structure, though they are present [2]. Therefore, for 
many problems and for many practical purposes, the folding configuration may be considered 
as topologically planar, which graphically amounts to the rule to draw the sequence and 
the pairings on the plane without self-intersection (figure IT]). This approximation makes 



* Electronic address: 



fr axLcois. david@cea. fr , wiese@lpt .ens. fr 



2 



the problem of RNA folding considerably simpler, since it allows for instance a recursive 
calculation of the partition function of a RNA strand in a polynomial time (as a function 
of the length of the strand). A lot of work has now been invested to find the most efficient 
algorithm [21 IH IS] • The planar approximation is also the starting point for the study of more 
general configurations, by performing expansions in terms of the topological number of the 
latter. Such studies may involve beautiful mathematical tools like random-matrix theory 
i6j. 

The folding of planar configurations of RNA strands is a fascinating subject in itself, 
with a lot of attention from physicists and mathematicians (besides biophysicists and bio- 
chemists). In particular, planar folded configurations are topologically equivalent to tree-like 
configurations, and the statistics and combinatorics of trees is a vast subject of its own. The 
homopolymer problem (all bases identical) was already solved in 1968 by de Gennes [7]. 
In this simple case the pairing-probability V of two RNA-bases decays with the distance ^ 
between the two bases along the backbone according to the scaling law V{t) oc Irre- 
spective of the embedding in 3-dimensional space, the statistics of the configurations is that 
of so-called "generic trees", or equivalently of the mean-field approximation for branched 
polymers. For Real RNA molecules however, the optimal fold depends on the sequence. 
Most studies, in particular numerical ones, focus on the configuration space and on the 
statistics and dynamics of folding for specific (and biologically relevant) sequences [21 18]. 

Since the pioneering work of Bundschuh and Hwa [Q UHl EH 112], several authors have 
studied the statistical physics of RNA secondary structures for random sequences and ran- 
dom bond energy models [THIIII]. One motivation is to understand the relative role of general 
sequence disorder and of specific biological sequences in the behavior of long RNA strands, 
and whether some properties are generic irrespective of the details of the sequence. In ad- 
dition the physics of random RNA sequences is interesting in its own, as a highly nontrivial 
example of (seemingly ID) disordered systems, where ordering (due to attractive pairing 
interactions) and frustration (due to the sequence disorder and the topological constraint of 
planarity) coexist. A key feature of the above models is that there appears be a continuous 
freezing transition between a weak-disorder phase, at large scales indistinguishable from the 
homopolymer case, and a strong-disorder or glass phase with non-trivial scaling, and of pos- 
sible biological relevance since the conformation and properties of RNA depends on the se- 
quence disorder, i.e. on the primary structure. Not much is known about the transition, even 
from numerical work; indeed its localization is non-trivial [11]. Better studied numerically 
is the glass phase at strong disorder, or equivalently zero temperature [HI UHl [HI [121 US US] . 
However, the nature of the freezing transition and of the low-temperature phase are still 
poorly understood, and contradictory results are reported [15]. It is e.g. disputed whether 
replica-symmetry breaking exists in the latter [161 [13 IIH]- The glass phase appears in the 
solution of [HI Uni [m [12] for the partition function for n = 2 replicas (instead of n = 
relevant for the disordered system) and in numerical simulations [HI [IHl El [121 [131 [U] • One 
feature which seems to be robust is the pairing probability V{t) oc i~Psi^^<> with Pgiass ~ |j 
independent of the disorder, be it sequence disorder, or random-pairing energies [IHl 111] . 

To better interpret the numerics, finite-size effects have to be understood. A first step 
in this direction was the recent analytic solution of a simplified hierarchical model [19], 
corresponding to a broad distribution of pairing energies, and with a pairing exponent p = 
(\/T7-3)/2 = 0.5615.... 

Pulling a DNA molecule at both ends has become an important experimental technique, 
which may one day allow to identify the DNA sequence by its force-extension characteristics 
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|20j . For RNA, the problem is more complicated, since folded RNA is not a linear strand, 
thus the sequence of base-pair openings is not clear in advance. There is a rapidly increasing 
bibliography on the subject [211 [22] • Remarkably, RNA pulling gives one of the first direct 
tests [2S| of Jarzynski's equality [211 [2S]- Averaged quantities can more easily be estimated 
and measured, either for homopolymers or numerically for disordered sequences [13] . Efforts 
have been undertaken to include experimentally relevant details, as the elasticity of the free 
RNA strands [261 E?]. 

B. The field-theory approach 

This paper is devoted to a renormalization group study of the freezing transition and 
of the force-induced denaturation transition of RNA with random pairing energies. In [28] 
Lassig and Wiese (LW) pioneered a field theoretical approach for the freezing transition 
for this model. They proposed a continuum formulation for the perturbative weak-disorder 
expansion of random RNA. Its starting point (the free theory) is the homopolymer model. 
They analyzed the divergences of this expansion at first order in the disorder strength, and 
they showed their model to be renormalizable at first order in perturbation theory. Assuming 
scaling at the freezing transition, they showed that this transition can be described by an 
UV stable fixed point at finite disorder strength, and that the couphng (disorder strength) 
and the length of the RNA strand (number of bases) have to be renormalized at one-loop 
order. This allowed them to calculate the critical exponents (to be described later) for 
the freezing transition. Using a "locking argument" (see below), the scaling exponents for 
random RNA in the strong disorder phase were estimated, in good agreement with numerics 
[91II011II11I211T3]. 

It is important to understand if this approach defines a consistent theory beyond first 
order (if possible to all orders), and if the estimates of for the scaling exponents are 
reliable. Indeed the diagrammatics in the LW model is of a new type, although it bears 
similarities with the diagrammatics of the Edwards model for polymers, i.e. self-avoiding 
random walks, and self-avoiding polymerized membranes, with non-local interactions. It is 
not at all obvious if the (now standard) field-theoretic renormalization formalism, leading 
to the renormalization-group picture, is valid for this kind of model. It is the purpose of 
this article to show this, and to present applications of this field-theoretical formalism. We 
shall introduce a formulation of the LW model in terms of interacting random walks in 
(i = 3 dimensions, and field-theory tools developed for self- avoiding membranes [221 EOl ED 
E21E31[SlE3[MlEZlEHlESllinillIlll21ll3]. We show that this model is consistent and 
renormalizable to all orders of the weak-disorder perturbative expansion, and deduce that 
the LW model is indeed renormalizable. Our formulation is in fact more convenient for 
explicit calculations than the original LW formulation. It allows us to derive new scaling 
relations between exponents, and to calculate critical exponents at second order. A short 
summary of this approach and its results at second order has already been published 
Our formulation allows also to treat the related problem of the denaturation transition of 
RNA strands induced by an external pulling force. The modelisation of this effect and the 
principle of the renormalization group calculation has been presented in [13] by the two 
authors and C. Hagendorf, but the details of the second order calculation are presented for 
the first time here. 
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C. Organization of the article 



The article is organized as follows: In section [TT| we discuss basic properties of RNA 
molecules and their folding, the equivalent description of these foldings in terms of trees, 



arch systems, and random-height models in subsection 11 A We then present the Lassig- 



Wiese field theory for RNA folding, firstly for the free theory (no disorder) in subsection 
11 B, and secondly for the interacting theory (with disorder) in subsection 11 C The per- 
turbative expansion of the interacting theory, its short-distance (UV) singularities and its 



renormalization are briefly discussed in subsection 11 D 



In section III we introduce our representation of the model in terms of interacting random 
walks. The basic idea relating random planar foldings to random walks in 3-dimensional 
space is recalled in subsection 111 A The representation of a free folded RNA strand (no 



disorder) in terms of a closed random walk and the precise concepts and notations are given 



in subsection 111 B We then generalize this representation to open random walks, since this 
will prove convenient for renormalization. In order to take into account the planarity of 
the folding, we introduce auxiliary "dressing fields", before taking a large- limit {N is 



the number of components of these fields), as detailed in subsection lllC The disordered 
( "interacting" 



model is introduced in subsection HID 



Since the disorder in the random 
pairing energies is quenched, we introduce n replicas. The average over the disorder gives 
an effective non-local interaction between replicas, given by the so called replica-overlap 
operator \t'. Finally, one takes the n — ^ limit ("replica trick"). The principles of the 
perturbative expansion and its diagrammatics are given. The model and its diagrammatics 



are easily extendible to an interacting open RW (open strand), subsection lllE, and to 



multiple interacting RWs (multiple strands), subsection lllF This is required in order to 
extract all renormalizations without going to 3-loop order. 

Section IV deals with the short-distance UV divergences, and their renormalization. 
Defining the model in (fictitious) d dimensions, with d = 3 relevant for RNA folding, dimen- 
sional analysis shows that d may be used as an analytic regularisation parameter (dimen- 
sional regularization) , and that the model is expected to be renormalizable for d = 2. This is 



briefly explained in subsection IV A The crucial tool to analyze the short-distance singulari- 
ties is the Multilocal Operator Product Expansion (MOPE), which generalizes the standard 
Wilson OPE for local field theories. It is an extension of the MOPE introduced in [311 132] for 
self-avoiding manifold models, and used in[3ll|32l|33l[3l[35l|36l[371[3Hl[33lig 
interacting tethered membranes and polymers. This MOPE and its structure for the differ- 



ent relevant (local and multilocal) operators is introduced and discussed in subsection IV B 



In subsection |IV C| we use the MOPE formalism to analyze the UV divergences of our model 
for random RNA folding, and show that it is indeed renormalizable (for e = — 2 — > 0). In 



subsection IV D we discuss the general structure of the counterterms and of the renormalized 
action. We show that UV finiteness requires a renormalization of the coupling constant g 
(as expected), a renormalization of the field r (which represents the position of the random 
walk in the ci-dimensional fictitious space), plus an additional renormalization for a bound- 
ary operator in the case of open RWs, which is crucial for the consistency of the model. The 
definition of the renormalization-group beta-functions and of the anomalous dimensions of 



the operators is given in subsection IV E Two slightly different renormalization schemes, 
denoted MS and MS, and based on the standard minimal subtraction scheme (subtraction 
of poles in e = d — 2) are introduced in subsection IV F They will be used for the explicit 



calculations. In subsection IV G we discuss renormalization of the so-called contact operator 
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and show that its anomalous dimension is not independent of the renormahzation of r, 
thanks to a new scahng relation that we derive with the help of the multi-strand model. 
Finally, in subsection IV H we apply our results to show that the Lassig-Wiese model for 
random RNA folding is indeed renormalizable, and we make precise the relation between 
our renormahzation of g and of r and the renormahzation of the coupling constant g and 
of the RNA strand length L in the LW model. This was the first initial motivation of our 
study. 

We then compute at second order (two loops) the renormalization-group functions of 
the random RNA folding model, and the scaling exponents for the freezing transition. In 
section |V] we give explicitly all diagrams and integrals. In subsection V A we present the 
principle of our calculation for RNA strands of fixed length. In subsection VB we present the 
calculation for another ensemble, the so called "grand canonical" ensemble, where the length 
of the RNA strand is a fluctuating variable distributed with an exponential distribution 
involving a chemical potential r. This is reminiscent of the two ensembles present for self- 
avoiding polymers (fixed length), and the n = field theory the after deGennes mapping 
(Laplace transform) ^6]. The calculations require the evaluation of diagrams to two loops 
in perturbation theory. This is done in subsection V C for one-strand configurations, in 
subsection VP for two-strand configurations, and in subsection VE for diagrams involved 
in the renormahzation of the contact operator $. 

In section VI we apply our two-loop calculations to the freezing transition. In subsection 
VI A we sum the results of section W\ and compute the ultraviolet poles in e = d — 2 for 



the partition functions at two loops, for an arbitrary number of replicas n. This determines 
the counterterms, the beta function and the a nomal ous dimensions, in the MS scheme (sub- 



section VI B), in the MS scheme (subsection VIC), and in the grand-canonical ensemble 
(subsection VI D). In subsection VIE we study the RG flow for e > 0. We show that the 
two-loop calculation confirms the existence of an UV stable fixed point (i.e. a phase transi- 
tion) at positive coupling, for n = as well as for n > 0. The n = fixed point describes 
the freezing transition induced by strong enough disorder in the random RNA model. We 
compute the critical exponents at second order in e, and check that the results are consistent 
between the different schemes and the different ensembles. 

Finally, in section VII| we generalize our approach to an applied external force pulling on 
the RNA strand. This problem was first studied by the two authors and C. Hagendorf in |15] . 
where it was shown that our model could be extended to describe the denaturation transition 
induced by an external pulling force, and where a one-loop calculation was performed. In 
subsection 



VII A 



subsection 



VII C 



we recall the model, and in subsection VII B its diagrammatics, while in 
we derive its renormalizability, and define the form of the renormalized 



action and of the RG functions. In subsection VII D we present new results, namely the 
details of the two-loop calculation of the counterterms and of the RG functions. In subsection 



VII E we discuss the physical meaning of our calculations for the influence of an applied force 



on the freezing transition, and on the nature of the denaturation transition for weak and 
strong disorder. 

Section VIII offers conclusions and further perspectives. 
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II. THE LASSIG-WIESE FIELD THEORY 
A. RNA folding representations 

1. Pairing configurations 

As explained above, RNA molecules consist of 4 bases, adenin, guanin, cystein and 
thymin, which are attached to a sugar phosphate backbone. In contrast to DNA molecules, 
there does not exist a complementary strand, and the RNA molecule has to fold back onto 
itself. For the reasons alluded above, we consider that an allowed RNA-fold is an RNA 
configuration, which can be drawn in the plane without self-intersections, see figure [l] (a). 
Equivalently, this can be redrawn as a set of arch- diagrams, given on figure [l] (b), or as a 
height diagram with the constraint that the height is zero at the both ends, and changes by 
zero or one between neighbors, figure^ (c). This will be explained below. 



(b) 



(c) 




s 



FIG. 1: (a) an (open) RNA molecule. Bases are represented as discs (color-coded for the bases); 
pairings are indicated with dashed black lines, (b) the same configuration, redrawn s.t. bases lie 
on a line; allowed configurations consist of non-intersecting rainbow configurations, (c) the same 
configuration, redrawn in the height picture: each time a circle in (b) starts to the right, h{s) 
is increased. Each time it comes from the left, h{s) is decreased. The such constructed height 
function h{s) of base s is the minimal number of bonds that have to be opened in (a) to reach base 
s. It is also equivalent to the number of arcs over base s in (b). 

In order to describe the LW model, we need to more precise. We consider a strand with 
L bases, i.e. of length L. We label successive bases by integers i = 1, . . . , L, and denote 
the pairing between two different base i and j by the ordered pair A planar pairing 
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configuration $ is given by the collection of N pairings 

* = {{il, k), {is, U),--- (^27V-l, i2N)} 

such that all the iaS are different and such that the corresponding configuration is planar, 
i.e. no knot or pseudoknot configurations are allowed. This implies that for any two pairings 
in # we have 

'^2 < is, 

ii < il, 

ii<is<i4< l2, 
<ii<i2< ii- 

A pairing configuration is compact if all bases are paired, that is if L = 2A'". 



1^1, , 1^3, «4j 



either 



(1) 



* * * 



II I i\i \ 

II I I 1^ \ 

i I ^ - / / M * ^ 



FIG. 2: Planar (left) and non-planar (right) pairings 



Any planar pairing configuration can be represented by an arch system. Associating to 
each interval i =\i,i + 1[ the number of arches i.e. the height h{i) above it, each planar 
configuration is in one-to-one correspondence with a path i — > h{\) over the non-negative 
integers with increment h{\. + 1) — h{i) = 0, ±1 for general configurations ("Motzkin paths") 
or h{i + 1) — h{i) = ±1 for compact configurations ("Dyck paths"). 

Finally, to each pairing configuration $ we associate the pairing function ^{i,j) which 
is defined by 

f 1 if (i, j) e ^, i.e. if i and j are paired, 
' 1 otherwise. 

Defining by r]{i,j) the pairing energy between to bases i and j, the energy of a folded 
configuration is 

m;v]= Yl (3) 

l<i<j<L 



2. Scaling exponents 

Irrespective of the precise statistics of pairings, these representations allow to define two 
scaling exponents, ( and p, which play an important role in the study of RNA folding. First, 
the average height (h) scales with the size L of the RNA-molecule as 

Jh) ^L^ . (4) 

Here, we denote by (. . .) thermal averages, and by an overbar . . . disorder averages. 
The probabihty that bases i and j are paired, scales hke 

WiJ))--\i-j\-' ■ (5) 
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provided that 1 ^ — j| ^ L. These exponents are not independent [28]. Note that the 
height at position k is the number of rainbow-arches starting before and ending after k 



(6) 



0<i<A: k<j<L 



Summing over all k on both sides and taking thermal and disorder averages, yields from 
scahng, assuming that 1 < p < 2 

L^+i ^ L^-p . (7) 
This yields the important scaling relation 



P + C 



(8) 



The pairing statistics depends on the set of pairing energies ri{i,j). For "homo-polymers", 
i.e. a uniform rj{i,j) = t] < for all i,j, de Gennes [IB] has shown that p = po = | and 
C = Co = ^. This can be understood from the fact that the height h{s) is a random walk 
in time s, constrained to remain positive (see subsection II Bl). These exponents are also 
relevant for random RNA in the high-temperature phase. 



We define the "pair contact" probability and the exponent 9 by 



(9) 



We expect that in the high-temperature phase = 2p = 3, whereas in the low-temperature 
phase this relation is not satisfied. In the glass phase, and if the partition function is 
dominated by a single or a few configurations, 6'giass = Pgiass- Since 



(10) 



it follows that in all cases 

p<9<2p . (11) 

We expect that upon lowering the temperature, there will be a phase transition with different 
universal exponents p* and 6*. Finally, in the low-temperature phase, there is a third set of 



exponents Pgiass and Cgiass- AH these exponents must satisfy relation (11). 



3. Free energy, finite-size scaling, and divergence of specific heat 

General scaling analysis yields that close to a fixed point g* of the renormalization-group 
beta function P{g), i.e. close to a phase transition 

Ig-g^l-r^'^'*^ , (12) 

where ^ is the correlation length. As will be explained later, the coupling g comes with the 
pair-contact operator, and thus 

P'ig,) ^2-9* . (13) 
Since close to the transition g varies continuously with T, this gives 
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The free energy scales like the inverse correlation volume, i.e. in one dimension like 1/^ ^ 

~ Lli . (15) 



Using (14), the divergence of the specific heat becomes: 

d^ , , 29*-3 



dT2 



J= ^ \T -T^\^ . (16) 



Thus for 1 < 0* < 3/2, this phase transition is of second order. 

For our model, it is difficult to extract the correlation length ^ from a simulation or exper- 
iment, since there is no scale at which a correlation function starts to fall off exponentially. 
Rather, ^ is the scale, where the contact probability ($(i,j)) ~ \i — crosses over from 
p = I to p = p* or p = Pgiass, which will turn out to be p* ~ Pgiass ~ |- 

B. The free theory 

We now recall the formulation of the LW continuum theory in the case where there is no 
disorder. 



1. Counting configurations 

In the absence of disorder, all pairing configurations are assumed to be equiprobable, with 
the topological constraint that they must be planar configurations. If we restrict ourselves 
to the case of compact configurations, the number of planar pairings for a strand of length 
L = 2N (number of Dyck paths of length 2A^) is given by the Catalan number 

Cn:= y l = ^(^^] ^ iV-^/^4^7^-l/^ (17) 

^ N +1 \ N J AT-^oo ^ ^ 

planar compact # 

In the general case (Motzkin paths), or in more realistic models where there is a weight for 
forming an arch, the number of planar configurations obeys a similar asymptotics 

Zo(L) = V 1= ~ L-P°c^ao with po = 3/2 , (18) 

' L— >oo 
planar # 

where Cq and Oq are non-universal constants. The exponent po = 3/2, which governs the 
power-law correction factor is a universal scaling exponent playing an essential role in 

the problem. The value of this exponent can also be understood from the observation, that 
in the height formulation of the problem, the planar pairing ensemble becomes a random 
walk ensemble on the half line h{i) G N. The exponent po = 3/2 is then nothing but the 
exponent for the probability for the first return to the origin at large time for a random 
walker in one dimension. 

^ A simple model for the last equation is as follows: Suppose that the system is correlated over a size 
^. Then there are {L/^Y independent uncorrelated degrees of freedom. If they have Ising character (2 

states), then (5 J" = -ksTlnZ ^ -/csTln ('2^^ = -ksT (^^ \n2. 
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Zo{L) is the partition function for planar pairings of a RNA strand with length L, when 
the energy for every possible pairing is the same, and when the only constraint comes 
from the planarity condition. This problem is similar to the problem of folding for an 
homopolymer considered by de Gennes in 1968 [7]. 

Let us now consider a strand with length L and impose the constraint that there are P 
fixed planar sub-pairings in the configuration. This collection of sub-pairings is denoted 

f^ = {(^l,Jl),(^2,j2),-- - ,{^P,JP)} ■ (19) 

These P sub-pairings divide the structure (with length L) into P + 1 substructures with 
backbone lengths Lq, Li, . . ., Lp such that 

Lo + Li + --- + Lp = L-2P . (20) 

The number of planar configurations with the substructure f2 fixed is denoted Zq{L\Q), and 
is given by the product of the number of configurations in each substructure, 

Zo(L|f2)=^l = ^<l>(^i,Z2)---$(zp,jp) = Zo(Lo)Zo(Li)---Zo(Lp) . (21) 

We write the sum over configurations as an unnormalized expectation value of an observable 
O as 

E O = (O)o ' (22) 

so that the partition function is the e.v. of the "unity operator" 1 

Z,{L) = {1), , (23) 

while the number of configurations with a fixed planar substructure f2 is the e.v. of the 
operator 

^{iu--- ,Jp) = $(ii,«2)---$(^p,Jp) (24) 

and reads 

Zo(L|0) = {Q), = mti,n)---<^{tp,jp)), (25) 

With these notations the number of configurations with a fixed planar substructure behaves 
in the large-size limit {L — > oo, all Ln/L fixed and of 0{1)) 

Zo{L\n) = ( ^^ )o ^ ^ L.^'L-^' ■ ■ ■ L-/^ c^^^ (26) 
2. Continuum theory 

In the rest of this article, we are interested in the scaling behavior for long RNA strands, 
and take the limit L ^ oo. We consider RNA in presence of disorder induced by the 
heterogeneity of the base sequence (primary structure), and will construct a perturbation 
theory in the strength of the disorder. In this perturbation theory, each term involves 
the expectation value for a product of a finite number of $(z,j). Our starting point is a 
continuum free theory where: 
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FIG. 3: Two planar configurations with the same fixed substructure fl, and graphical representation 
of Z{L,ft). 



1. The length of the strand L is rescaled to be finite. 

2. The positions ia of the bases become a continuous variables Sa, 

2G{1,2,---,L} ^se[0,L]. (27) 

3. The non-universal factor CqOo in the partition function Zq{L) is absorbed in the nor- 
malization of the expectation value ( )q so that the continuum partition function of 
the strand with length L is 

Zo(L) = ( 1 )o = , po = 3/2 (28) 

4. Similarly, the non-universal factor ao/cg is absorbed in the normalization of the oper- 
ator $(i,j), s.t. in the continuum limit the operator fl is defined by its expectation 
value 

Zo{L\n) = {n), = ms^,t^)■■■^sp,tp)), (29) 

with 

, ^, \ Ln^" L7^" ■ ■ ■ L'p^" if f2 is a planar structure , 

($(si,ti)---<l>(sp,tp))o = ' . (30) 

I otherwise. 

3. Diagrammatic representation 

A convenient diagrammatic representation of the operator $ in the free theory is the 
following. We represent the partition function for a strand of length L, Zo{L), by a single 
line with length L. This single line represents the whole (normalized) sum over all planar 
pairings between points on the line. 




(l)o = = E ^ • (31) 
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The operator $(s, t) is then represented by a dashed arch over the hne joining points s and 
t (it is a bi- local vertex joining s and t): 



{Hs,t)),= ^ ~ . (32) 

The partition function for a strand with fixed planar substructure f2, Zo{L\n) = {^1)q, is the 
expectation value of a product of $ operators and is depicted by the corresponding planar 
collection of arches over the line. If the substructure Q is non-planar, the expectation value 



(i7)o is zero, according to (30). 




{n)o = 1 ' = L ^ ^ . (33) 

<s>dq 

The expectation value of fl depends only on the sub-backbone lengths, hence on the distances 
between the end-points of the arches considered to be on a closed circle. Both endpoints 
of the strand are identified, since this does not change the statistics. There is formally no 
difference between open and closed RNA strands, since we are interested in the secondary 
structure, not in the tertiary structure, thus steric effects are absent. An alternative dia- 
grammatic representation for the partition function and the $ operators is to depict Zq as 
a closed loop with a marked point which depicts the endpoints of the strand. Similarly, 
the partition function for a strand with a fixed planar substructure Q, Zo(L|f2) = (fi)^ is 
depicted as a closed planar arch system. This is represented on figure |4} 



C. Random RNA, disorder and the Lassig-Wiese field theory 

1. Random RNA 

The field theory approach initiated in [28j by Lassig and Wiese is based on the random 
RNA model proposed by Bundschuh and Hwa [HI E]. In this model one assumes that 
to each pair is associated a pairing energy T]{i,j), and that the total energy E for a 
pairing configuration $ is the sum of the pairing energies associated to each pair. With our 
notations the configurational energy E may be written as 

E[^-r]]= v{^,j)Hhj), (34) 

^<i<j<L 



with $(i,j) the contact function defined by (43). Given the collection of pairing energies 
V = {vihj)] ^ < j}) the partition function for the RNA strand at finite temperature is the 
sum over all planar configurations 

Z, = ^exp(-/5E[$;r/]) , (35) 



13 




FIG. 4: Equivalent diagrammatic representation as closed arch structures. 



with 13 = l/ksT the usual Boltzmann factor. For fixed pairing energies t], the partition 
function of a strand of length L can be computed recursively (see e.g. [HI US]) in a time 
T = 

While biological sequences are highly structured in order to fulfill their biological function, 
here we consider random RNA sequences. While this may, or may not be realistic for real 
RNA, it is at least an important benchmark against which to compare experimental results 
for biologically functional RNA. 

However, the random-sequence model is still not amenable to an analytical treatment. 
We therefore assume that the pairing energies are independent Gaussian random variables. 
This approximation neglects correlations between the random pairing energies. Numerically 
it seems that at least in the low-temperature phase, these correlations do not affect the 
large- distance properties [13]. We have to leave to future research to develop an analytical 
handle on this problem. 

We choose to be a random variable with probability distribution 

r^o < is the mean pairing energy and a its variance. The randomness in the pairing-energy 
distribution amounts to the introduction of quenched disorder in the system. The average 
over the Gaussian disorder t] is denoted by the horizontal overline , so that (with the 
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ordering i < j and k < I) 

r]{ij)ri{kj) = 7]q + a'^ 6ik6ji . (37) 
The averaged free energy for the system is 

^ 1 



F = --logZ, , (38) 



and the expectation vahie for an observable O is 



(0):=^J]Oexp(-/?E[$;r/]) . (39) 



For instance the probabihty that the bases i and j are paired is 



($(2,j)) = ^$^<^>(^,j)exp(-/?E[$;r/]) , (40) 



and the probabihty for a given pairing substructure Q to occur in the random pairing-energy 
ensemble is 



(fi) = ^5^fiexp(-/3i?[$;r^]) = exp(-/3E[<|.;r/]) . (41) 



2. Weak-disorder expansion and replicas 



The idea of [2H] is to study the model by a perturbative weak-disorder expansion, and 
to extract its large-length L scaling behavior at finite (and if possible large) disorder by 
renormalization group techniques. The perturbation expansion can be constructed in the 



discrete model by expanding in powers of the disorder //(i, j) and using (37) so that we get 



a perturbative expansion in powers of the effective disorder strength (coupling constant) 

g = P'a\ (42) 

The quenched average over the disorder is done by the standard replica trick, which is well- 
defined for a perturbative expansion. One considers n replicas of the system, labeled by an 
index a = 1, 2, ■ ■ ■ , n. Finally, one has to take the limit ra — »• 0. The pairing configuration 
$a of the replica a is given by the pairing function $a(i, j) 

^ , , I 1 if (i, j) G $Q, i.e. if i and j are paired for the replica a, , , , 

Mhj) = \^ ' \ ^ ' . (43) 

I U otherwise. 

Since the disorder is quenched, all replicas see the same pairing energy ri{i,j), and the 
configurational energy for a replica ensemble $ = is 

n n 

S[<^;V] = Y.E[<I>^;V] = Y, E vihJ)<^a{^,J) • (44) 

a=l a=l l<i<j<L 
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The average over the disorder gives the partition function for the n-times rephcated system 



Z = {Z,r=Y, exp(-/3^[<|.;r/]) . (45) 

The average over the disorder can be taken exphcitly since the disorder is Gaussian. From 
(37) one has 



exp(-/35^5^r7(2,j)$a(^,j)) =exp(| Y.^'^^'^^^^P^'^^A ^^^^ 



Q=l i<j / \ a,/3=l i<j 



with g given by (42). One obtains an effective attractive interaction between rephcas, and 
one can rewrite the system in terms of an effective "Hamihonian" 

n n 

/37^tot[<^>]=/9^o[<f]+/9^[<f]=/3r7o$^X;$a(^,J)-f 5^ 5^$.(^,J>/3(^,J) • (47) 

a=l i<j a, 13=1 i<j 

The first contribution, proportional to (3t]o is the one present for a homopolymer, which 
we can solve analytically. The second term, proportional to the disorder, contains two 
contributions: the diagonal contribution a = P 

i<j i<j 

leading to a change of rjo 

n 

and an off-diagonal part 7Y[$] 



a=l i<j 



7^[$] = -^5^^^,^(^,j) , ^,^(^,j) :=<fa(^,j)'^'/3(^,j) , (50) 

a</3 i<j 

where \i/ 0,^(2, j) with a ^ (3 is the pair contact or overlap operator. It gives the probability 
that the bases i and j are paired both in replica a and (3. The partition function and the 
e.v. of observables are now 

<!> " " \ <!> / 

(51) 

The idea is to do perturbation theory in H, using the solvable theory with TCq as reference. 



3. Continuum limit 

The Lassig-Wiese field theory [28j is obtained by taking the continuum limit of this 
model. It is defined in terms of the continuum pairing operators for each replica ^a{u,v), 
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M, f G [0,L], and the overlap or pair-contact operator \I/q,^('U, f ). One starts from the free 
theory for n independent non-interacting rephcas. The partition function for a bundle of n 
free rephcas is 

Zo{L) = J2^ = Il^^)o = Zo{Lr = L-"'^^ . (52) 

We represent diagrammatically this partition function by a collection of n lines, or by a fat 
bundle. 



2o = = n -^^^=^^^=' (53) 

a 

The expectation value for a product Q of n (different) operators Qa living in replica a 
factorizes into 

Zo{L\n) = J2Il^c. = mo = ll{^a)o = llZo{L\^a) • (54) 

{$„} a a a 

We represent it diagrammatically by the collection of the n planar arch structures relative 
to each Qr, 



1- I 

I ^ L_ 



Zo{L\n) = ' ' ; ^ . (55) 

The continuum model with disorder is given by the theory with an effective disorder Hamil- 
tonian corresponding to an attractive 2-replica interaction which is the continuum limit of 



the discrete effective Hamiltonian (50) 



a</3 



nm = -gJ2 '^ap{u,v) , '$o.p{u,v) ■.= <^^{u,v)<^p{u,v) . (56) 



0<u<v<L 



The partition function for n replica of a strand with length L with disorder is the continuum 



version of ( 45 ) . Therefore it is given by 

ZiL) = {expi-n[mo ■ (57) 

It will be expanded in powers of the coupling constant g, each term of order g'^ being of 
the form {n[<l>] )q and can be computed using ( [54| ). Similarly the partition function for n 
replicas with a given set of substructures Q (i.e. the "expectation value" for the operator Q) 

with disorder is 

Z{L\Q) =l^r) = {Q exp{-nm))Q . (58) 

It will be computed as a formal power-series expansion in g in terms of the )o. The 

details of these perturbative expansions will be discussed and studied in the next sections. 
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r s t uv wx y l 

FIG. 5: Example of a product of four ^'-operators, consistent with planarity. 

D. Perturbative expansion for the Lassig-Wiese theory 

According to the LW diagrammatics, the continuum overlap operator \l/Q,^(M,f) is repre- 
sented by a double arch between points u and v on the two lines for replicas a and (3 

^o.f3{u,v)= / ^\\ . (59) 



The perturbative expansion in g of the partition function 2{L) involves integrals of e.v. of 
products of \E' operators 

J J J {Ui,Vi} 

which is represented as a set of k double arches between the replicas. An example is given on 
figure [s} At a given order g'^, the number of different replicas coupled by the \l/-arches ridis is 
bounded by 2 < < 2A;. Let us consider a configuration \[' ■ ■ ■ \[f associated to replica pairs 
{ai,Pi) and base pairs {ui,Vi); for each replica among the ridm coupled replica, and consider 
the corresponding reduced system of $-arches. The e.v. (\E'q,^^^(mi, fi) ■ ■ ■ 'ifakiSki'^kyVk))^ is 
non-zero if and only if for each replica a the reduced system of ■ ■ ■ is planar. For 
each replica a the e.v. of the product $q, ■ ■ ■ is the product over the + 1 cycles of their 
backbone lengths, to the power — po = — 3/2 

($„---$,)o= I 174 J (61) 




The e.v. of the product ■ ■ ■ \[' is now the product of the previous terms for each of the 
ndiff coupled replicas, times the product of the free partition function = L~^/'^ for the 
n — ^difi uncoupled replicas. 

"diff / Pa \ ~p'^ 

(^:^>o=n n^^- (L--)"- (62) 

k a=l \j=0 J 



An example is given on figure |5] 
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We expect the scaling dimensions of $ and \1/ to change in the presence of disorder. 
It is the aim of this article to calculate these changes. In a perturbative field theory, the 
latter can usually be extracted from the divergences of the diagrams, as e.g. the one given 
on figure [5j We therefore have to achieve two things: Calculate these divergences, but 
even more importantly, find which quantities they renormalize. In a standard field theory, 
this task is not difficult: The needed renormalizations are associated to the marginal and 
relevant operators present in the original theory, or generated by the perturbation expansion. 
Here, and up to now, we only have a perturbation theory, but no field theoretic action to 
renormalize, so we do not know which quantities will need renormalization! In the following, 
we will construct such a field theoretic representation, which will tell us which quantities to 
renormalize. In a second step, we will then calculate the necessary diagrams. 



Before doing so, let us as an example calculate the diagram drawn in equation 59 , to see 
that indeed there are divergences: 



— Jo A ' ' ^ ' 2 r(2-4po) 

U V 

(63) 

The diagram has a pole in 1/(1 — 2po), renormalizing the free energy. It has also a pole 
in 1/(1 — po); the latter can be interpreted as a renormalization of the length of the RNA 
molecule. However, it is not at all obvious why, and how to do this, thus a proper represen- 
tation as an action is necessary. 



III. THE RANDOM- WALK REPRESENTATION 
A. Basic ideas 

The diagrammatic expansion of the model bears strong similarities with the diagrammatic 
expansion of the Edwards model [IZlllH], which describes 3-dimensional random walks with 
a weak repulsive interaction upon contact, and which has been widely used for polymers 
and self-avoiding membranes. 

Here is an heuristic explanation for this similarity: Planar pairing configurations for a 
RNA strand are in one-to-one correspondence with planar arch systems over a linear strand, 
which are themselves in one-to-one correspondence with discrete paths (Dyck or Motzkin 
paths) on the half-line of integers N. In particular the (normalized) partition function of 
the free strand with length L is nothing but the probability of first return to the origin at 
time t = L for a random walk on N, or Z, which scales at large times as the continuous 
random walk on M (the Wiener process), that is as Now for several observables, the 

one-dimensional random walk on the half-line IR+ behaves as the three-dimensional random 
walk on the full space M^. In particular, the first-return probability to the origin for a RW in 
one dimension scales as the total return probability to the origin for a RW in three dimensions. 
As a consequence, the pairing operator $(z, j) has a natural representation in the 3d RW 
picture as the so-called contact operator 6^{r{u) — r{v)) (probability of contact at times u 
and V for the random walk r(t)). Similarly, many observables and many questions about 
scaling can easily be represented or formulated in the RW picture. In particular, the analytic 
regularization used in |28], where the contact exponent po = 3/2 is analytically continued 
to p G [1,3/2] and used as an UV regularization parameter to construct an e-expansion 
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for the RG equations and the scahng exponents, is nothing but the classical dimensional 
regularization where the dimension = 3 of space for the RW is analytically continued to 
d G [2, 3], and used to construct a. d = 2 + e expansion. 

There is however an important difference. The planarity constraint for the pairings im- 
plies that the product of several pairing operators vanishes if the resulting configuration is 
not planar. This is a global topological constraint that cannot be represented by local oper- 
ators in a RW representation. To implement this constraint, we shall introduce additional 
matrix-like degrees of freedom in the RW representation which allow to deal with the topol- 
ogy of the diagrammatics and to take the planar limit as a large- limit (where N is the 
dimension of the "internal space" associated to these additional degrees of freedom). This 
is a usual trick in QFT and in statistical mechanics. In particular it has been introduced in 
^ for the problem of RNA secondary structure enumeration and statistics. 

Thus we construct in this section a quite involved RW-like representation of the LW 
model, which involves "generalized random walks" in a d x N x n dimensional space, where 
c? = 2/9 is the dimension of space, N is the dimension of internal space, and n the dimension 
of replica space. We are interested in the limit c? — >■ 3 (po = 3/2), N oo (planar limit) 
and n — > (limit of quenched disorder). 

This representation turns out to be very powerful. It allows to apply the mathematical 
tools developed in the renormalization of polymers and self-avoiding membranes, in partic- 
ular the so-called Multilocal Operator Product Expansion (MOPE). Although the random 
RNA model is mapped only on the closed RW subsector of the RW model, there are other 
observables, associated to open random walks, which have no interpretation in terms of RNA 
observables, but which are much easier to study and to compute. They allow a more direct 
calculation of some of the renormalisation group functions and of the scaling exponents for 
the random RNA model. 

B. The simple RW model 

In the rest of this article, we normalize the Dirac "5-function" in as 

S\r) = (47r)^/' 5''(r) . (64) 

With this normalization most of the annoying factors involving powers of An disappear in 
the calculations. 



1. Closed Random Walk 

We start from the random-walk process in M'' in the time interval t G [0, L], described by 
the random variable r{t) (r = {r^; n = 1 ■ ■ ■ d}). The Euclidean action for the RW is (with 
proper normalization) 

5o[r]= fdtU{t)\ (65) 
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and the functional measure is the standard Feynman-Kac measure ^ (in Euchdean time) for 
the quantum particle with mass m = 1/2 

= n . « = UV regulator ^ ^ . (66) 

The partition function for the closed RW (periodic boundary conditions) is defined as 

^(ciced)^^^ = I V[r] exp(-5o[r]) 5'^(r(0) - r(L)) . (67) 

To extract the infinite factor from the translational zero mode in M'^, we formally write 

Vol(M'^) = y"d'^r = (27r)'^5'^(q = 0) (68) 
(with q the momentum in the conjugate space of M"^), and 

^(closed) ^ ^^^y ^d^^ ^ ^(clccd) _ ^gg^l 

The normalized closed partition function Zq^^°^^'^^ (L) is nothing but the heat kernel given by 
^(closed) = (0| exp (LAr) |0) = / V[r] exp(-5o[r]) ~6\r{0) - r(L)) = L-"/' . (70) 

Jr(0)=O 



Note the disappearance of the usual (Att) '^I'^ factor, thanks to the normalization (64) for 
the Dirac 8 distribution in the definition of the partition function. 



For (i = 3, Zq'^'°'^'^*^'' (L) is equal to the partition function Zq{V) of the free closed RNA 



strand in the continuum limit of the LW model, as defined by (28) (hence the similar 



notation). Therefore we also denote it as the unnormalized expectation value 

^(closed) _ ^^^closed ^ (7^) 

and represent it as a single line with length L, as for the RNA model. 



^(closed) = = rT^ C^C:;-! (72) 




The normalization for the action (65) was chosen such that the propagator (the IR-finite 
2-point function) becomes 

l\r^[u)-r^{v)\Y{u)-r\v)\)^ = b^''\u-v\ . (73) 



The measure has dimension [i]'^/^ ~ [r]'^. 
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2. Open Random Walk 

Although the observables for the random RNA pairing model are related to observables 
for a closed random walk, we also consider open random walks with free boundary conditions. 
For the open RW it is convenient to define the generating function 

^^'"^(qi, q2; L)^ J V[v] exp {-So[v]) exp (i(qir(0) + q2r(L))) . (74) 

It is the Fourier transform of the partition function for an open RW with fixed boundaries, 

4°P^"Vi, L)^ [ V[v] exp (-5o[r]) . (75) 

Jr(0)=ri,r(L)=r2 

qi and q2 are the momenta flowing through both end points t = and t — L of the RW. 
With our normalization for the measure, and using translational invariance, it is given by 

4°P^")(ri,r2;L) = (ri|exp(LAr)|r2) 

= (47rL)-^/^ exp (-^^^^^) (76) 

hence _ 

^"'''^^(qi, q2; L) = (2^)' 5'{q, + q,) Zt'"'\cl^■, L) (77) 

with _ 

4°''"^(qi;i^)=exp(-Lq?). (78) 

These notations will be useful later. We represent diagrammatically the open RW function 
as a single line with length L with bars at its end points (if necessary for clarity). 




3. Contact operator 
The contact operator S(M,f) is defined as 

S(k, v) = ¥{r{u) - r{v)) = {ATrf' 5\y{u) - r{v)) . (80) 

Again, the factor of (477)*^/^ is a normalization factor simplifying the calculations. We repre- 
sent it diagrammatically as an arch joining the points u and v. The partition function with 
one contact operator inserted is thus 

/"/ \\closed ^, 

{^M), = 

= / V[r] exp(-5o[r]) 5(m, v)^{L-\u- v\)-'^'^ \u - v\-'^'^ (81) 

ir(0)=r(L)=O 
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This is nothing but the product of the sizes of all loops, raised to the power of —d/2. More 
generally, the partition function with K contact operators inserted is 

(S(Mi,fi) ■ ■ •S(M/^,fx))o°""^ = / e^v{-SQ[v])'^{ui,Vi) ■ ■ ■'"{uk.vk) 

Jr{0)=r{L)=O 

= PlKv,]-"^' (82) 

where Pl is the Symanzik polynomial of the K + 1 loop (0^-like) diagram obtained by 
contracting to a 4-vertex each arch associated to a contact operator, or equivalently the 
product of all loop sizes, raised to the power of —d/2. An example of such an arch system 
(with K = 10) is depicted in figure M together with the corresponding (f)^ diagram. 




Ui U2 U4 f 4 f 2 ^5 Uq Vq Vi Uj Ug Ug Vg Vg UiQ f lo Vl L 




FIG. 6: An arch system (top) , and the corresponding (f)^ cactus-hke diagram (bottom) with its 
loops. The length of the k = 11 loops are Iq = L — (vi — ui) — {v-j — uj), li = {u2 — ui) + (ns — V2) + 

(fl-Us), ^2 = {u-i-U2) + {Ui-V:i) + {v2-VA), £3 = {V3-U3), li = {Vi-Ui), 4 = {uq-u^) + {v5-vq), 
4 = {vq - Uq), £7 = {U8 - Uj) + {Uio - Vs) + {V7 - Viq), £s = (ug - Ug) + (^8 " Vd), ^9 = {V9 " Ug) 

and £io = (fio — uiq). The Symanzik polynomial Pl is the product of the length for each loop 
Pl = ni=o,io^* 



The main problem with this simple RW representation it that the product of several 
contact operators A does not vanish if the corresponding arch configurations are non-planar. 
As a consequence, the model of interacting RWs with a contact interaction given by the 
action 

S[r] = \ [ i{tf-9 II Hr{s)-r{t)) (83) 

^ Jo<t<L J Jo<s<t<L 

(the attractive Edwards model) has a perturbative expansion which contains, besides the 
planar contributions which correspond to terms in the expansion for a RNA pairing model, 
many more non-planar contributions. Thus the complete action must be more complicated. 
(Note that (83) is of course only a toy model, since it does not contain the replica part 
necessary to treat the disorder.) 
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C. The dressed planar RW model 

In order to classify the pairing configurations according to their topology in the RW 
representation, and to keep only the planar configurations, we introduce additional matrix- 
like degrees of freedom and modify the action accordingly. 

1. Auxiliary fields 

First we add a conjugate pair of auxiliary A^-component fields 7a (i) and 7a (t) with a 
dynamical Ito like action 

N „ 

Ut)ia{t) (84) 



'5o[7,7] = J2 

a=l JO<t<L 



a = 1, . . . N is a color index which will play its role later. The action is such that the 
propagator for these auxiliary fields is the causal Heaviside function 9 

I (J otherwise , 

while 

{la{h)%{h))o = (7a(ti)76(i2))o = V a, 6, ti, t2 . (86) 
We represent the propagator (85) by a dashed oriented line, see figure [7j The (closed or 



open) partition function for the free dressed RW is now defined as 

ML) = J V[r] P[7, 7] e-5«W-'5o[7-7] (87) 
where we insert the bilocal boundary operator 

^J2'^a{0h^{L)] (88) 




with the proper boundary conditions for the path integral over r(t). Fb creates an auxiliary 
field at the initial point and annihilates it at the end point, thus still giving a contribution 
1 for the free RW. At that stage nothing changes for the expression of the closed and open 



RW partition functions and Zq^*'", which are still given by (69) to (71) and by (76) 



to ( 78 ) respectively. However, the boundary operator is crucial for the correlation functions 
and the interacting theory. 

The dressed RW partition functions are now graphically represented as a ribbon (or fat 
line) with a full line for the RW and a dashed oriented line for the auxiliary field, see figure [7| 



2. Dressed contact operator 

Now we can dress the contact operator S(u, v) with the auxiliary fields so that it becomes 
a ribbon arch with its topological features encoded by the auxiliary field color indices a e 
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mih{t2))o = h ^ t2 



^(closed) 



FIG. 7: The auxiliary field propagator and the dressed closed and open RW partition functions 



{1, ■ ■ ■ , A^}. Let us define the dressed contact operator $ as 

Hu^^) = J^Yl ^"(^) ^biu) 5'^(r(n) - r(^;)) j,{v) U^) • (89) 

a,b 

With the diagrammatic rules given above it is depicted as a ribbon arch (see figure |8]). 



///■ 

' / ' \ \ \ 

a I I L._.b .0 1 L._.fl . 



U V 
FIG. 8: Diagrammatic representation of the dressed contact operator $(n, f) 



The partition function with the insertion of operators $, defined as in (82) by 
($(«!, V,)--- <!>{UK, vk))"^"'"" = [ V[r] P[7, 7] e-^«H-^"[^'^] $ («i, ^1) ■ ■ ■ $ («x, vk) 

7r(0)=r(L)=0 

= N-^'^ Pl[u,,v,]-''/'' (90) 



now contains a multiphcative color-counting factor N^"^^ where h is the number of handles 
of the surface on which the planar arch system for the product of the $'s can be drawn 
without crossings. The argument is standard and requires counting the factors of which 
arise when summing over the internal color indices a carried by the auxiliary-field line, and 
by using the Euler relation for the system of arch diagrams. This counting is illustrated on 
figure [9] 

3. The planar N ^ 00 limit 

If we take the planar ^ 00 limit, only the contributions of the planar arch configu- 
rations (with h = handles) survive. This shows that we can build the RNA perturbation 
theory in terms of a self-avoiding polymer model embedded in = 3. In the following, 



we discuss how this can be put to work for the random RNA model defined in Eq. (47) 
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b • ' :' c \ 1 \ h 



X V 



u w 



V X 



FIG. 9: Example of planar (left, /i = 0) and non planar (right, h = 1) products of $'s. The left 
one carries an additional factor of A^^ from the sums over h and c. 



The advantage of this formulation, and the only reason we have gone through this formal 
exercise, is that we can write perturbation theory with a polymer-like action [microscopic 
free energy), which allows us to apply the tools of non-local field theory, and especially the 
multilocal operator product expansion MOPE (see [33] for a review). 



D. Replicas, the interacting RW model and its diagrammatics 

1. The action 

To construct a random walk representation of the LW field theory, we introduce replicas 



and construct a RW representation for the effective interaction term (56) between replicas 



induced by the quenched disorder. Consider n replicas of the RW field, rQ,(t), labeled by 
a G {1, ■ ■ ■ and n replicas for the auxiliary fields 7"(t) and 7"(t). The action for the 
free replica system is 



Jo<t<L ^ „, „ „ 



The dressed contact operator for the replica a is 

a,b 

with 'Ea{u,v) the contact operator for the replica a 

The dressed overlap operator between distinct replicas a ^ (3 is 

^'a/3(M,f) = ^c,{u,v)^p{u,v) . 

The attractive replica interaction is 



'5int[r,7,7] = -9^ 



0<u<v<L 



SO that the full action is 



S = Sq + S; 



(91) 



(92) 



"T '->int • 



(93) 



(94) 



(95) 



(96) 
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FIG. 10: Diagrammatic representation of the contact operator (left), and of the overlap operator 
^afs (right) 



Zo = = — 

FIG. 11: The replica bmidle. 

We shall construct the perturbative expansion and its diagrammatics for this theory, starting 
from the diagrammatic representation of the interaction operators represented in figure [TOj 
We are interested in the double limit — oo (planar diagrams) and n — >■ (quenched 
disorder). In the remaining discussion and in the calculations we shall take first the planar 
limit (the number of diagrams is thus greatly reduced), but keep the number n of replicas 
non-zero. We take the n — limit at the end of the calculation. Since = oo, we do not 
represent the auxiliary field propagators as dashed lines any more, but simply keep planar 
diagrams. 

2. Closed-strand partition function 
The n-replica closed-strand partition function for the free model {g = 0) is 

^(ciosed)^^^ = = [ P[r,7,7]e-^»rb = \ML)T = L-^"/' , (97) 

where the measure and the boundary operators are 

a a a \ a 



(98) 



We represent it diagrammatically as a single (fat) line, but it is understood that this repre- 
sents a bundle of n lines. 

The closed-strand partition function for the interacting model is 

^(closed) = (exp (-5i„t))S'°^^' = [ V[r, 7, 7] e"^ (99) 

^r„(0)=r<,(L)=O 

It can be expanded in a perturbation series in powers of g. The term of order g^ is of the 
form 

]^E///(*-«>o. (100) 
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and can be represented diagrammatically in terms of double arch systems involving P replicas 
with 2 < P < 2K, exactly as for the LW model. In the following we represent as a line 
only the replicas coupled via overlap operators. It is understood that the n — P other 

replica are there and give a factor of (Zq) " 

It is clear that for closed RWs the diagrammatics and the resulting integrals are equivalent 
to those of the LW theory. The amplitude ■ ■ ■ is given as a product over each of the P 
replicas a of products of the internal loop lengths to the power —d/2, equivalent to formula 



(I62j): 

P / Pa \ ^'^/^ 

(^)o=n ri^- (L-^/^)"-^ (101) 

k a=l \j=0 J 

The combinatoric factor for each amplitude, which is a polynomial of degree P in n, is 
also the same for the LW and the RW models. An example of such an amplitude is given 



in figure 12 The RW representation in the planar limit provides a (somewhat formal but 
systematic) functional integral representation of the LW model and a way to study its short 
distance structure (see next section). 



^ 'I 

" II 



r s t uv wx y L 

— [(^ ~ i^){y ~ — V + r — y + x)]^'^^'^ 
X [{u - t){v -u + t-r){L-v + r)]"''/^ 
X [{u — t){w — u + t — s){y — x){L — y + X — w + s)]~'^^'^ 
X [{w - s){L - w + s)]-'^/'^ 

X L^^-")"'/^ ^102) 



FIG. 12: Example of e.v. of A: = 4 ^''s as given in (101). The factors are grouped by replica lines, 
starting at the front. 



E. Single open-strand partition function 

As already explained above, it will be useful to consider other sectors and other observ- 
ables in the RW model. These observables are associated to open random walks and have 
no interpretation in terms of the random pairing RNA model. 

By similarity with the partition function for an open RW ( [74|[75 ), we consider now 



n 



replicas of a dressed open RW ra{t) with t G [0,L], and with free boundaries ra{0) and 
ra{t). We attach to each endpoint momenta qi and q2 which (for simplicity) are taken to 
be the same for the n different replica. The open-strand partition function is defined as 

Z(i)(qi,q2;L) = [[ P[r,7,7] e""^ Tb e^t^^^ '■"(°)+^^ (103) 



28 



where S is the action for the interacting dressed open RW model (96) and Fb the boundary 
operator, as defined by (98). The index added to the partition function (instead of 
the index (°p^^) used in Sect. IIIB IIIC) indicates: (i) that we deal with open strands, 
(ii) that we deal with one (1) single bundle of n replica of the same open RW. Later on 
we shall consider the partition functions Z^^^ ioi p > 1 (bundles of n replica of) open 
RWs interacting via the disorder-induced 2-replica contact-interaction term Using 
translational invariance in M*^ for each replica, we can factor out a momentum conservation 



term for each replica, thus defining, by analogy with the free case (77) 



Z«(qi,q2;L) = [(27r)'^5^(qi + q2)]" ZW(qi;L) . 



(104) 



The IR finite function 2'^'(qi; L) can be expanded in perturbation theory in a power series 



m g, as 



ZW(qi;L) = (e 



-iqi EJr„(0)-r„(L)]^°P«'i 



-5int^iqiEJr,(0)-r.(L)]\°P<''i 



(105) 



ra(0)=0 



with iSint given by (95). The rules to compute the perturbative expansion are a simple 



generalization of those for the closed RW model. The term of order can be expanded into 
a sum (over the various distributions of (a, /?)) of integrals (over the u and v) of expectation 
values of the form 



open 




(106) 



The integral over the auxiliary fields selects the arch configurations which are planar for 
each replica, and give zero for the others. We end up with a product for each replica of an 
e.v. for the open RW model of the form 



{A{ui,vi) ■ ■■A{up,Vp)e 



-iqir(L)\°Pe° 

/o 



(107) 



with the planar arch sub-system Q = A{ui, vi) ■ ■ ■ A{up, Vp) extracted from the planar dou- 
ble arch system '^aii3i{ui,vi) ■ ■ ■ '^aki3ki'^k,Vk) for each replica a. The e.v. (107) is easily 
calculated. The p arches form p internal loops with backbone lengths ii, - ■ ■ ,ip. At vari- 



ance with the closed RW model, the remaining segments of the strand that are not under 
an arch form an open sub-strand with total length io = if^-ee (with io + ii + ■ ■ ■ + ip = L) . 



The e.v. (107) is 



{A{ui,vi) 



A{up,Vp)e 



^iqir(L)\°P'=n 

/o 



n 



~d/2 



closed loops 



(108) 



and is represented in figure 13 



Each e.v. (105) is the product over P amplitudes of the form (107) (for the P replicas 



coupled by the \l/a/3), times the n 
by an amplitude of the form 



P remaining free-strand amplitudes. Hence it is given 



-cZ/2 



n n 



E io,c+in-P)L 



(109) 



a=l 



and is represented as in figure 14 
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{A{ui,vi) ■ ■■A{up,Vp)e 



-iqir(L)\°Pen 
/O 



V IV y 



Z X 



FIG. 13: Diagrammatic representation for the open strand planar arch sub-system (107) for a 
single replica, and the corresponding amplitude. The open strand with length £q is depicted by a 
bold line. The cut at the ends of the strand indicate that we deal with a open strand with incoming 
momenta qi. 



[{v — r){y — x)] X [{u — t){v — u + t — r)] 2 

_d _d 

X [{u - t){w - u + t - s){y - x)] 2 [(u; - s)] 2 

FIG. 14: Diagrammatic representation for a single planar open-strand ^ system, and the corre- 
sponding amplitude. The bars at the endpoints of the bundle indicate that we deal with a open 
strand with incoming momenta qi and — qi. 
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F. Multiple open-strand partition function 



Finally we consider partition functions for M > 1 strands interacting via the disorder- 
induced contact operator \E'. Let us restrict ourselves to the two-strand case M = 2. These 
two open strands are described by the RW's ri and r2, and more precisely by 



ria(ti), 7?Jti), 7i"Jti) (strand 1) 
r2a(t2), 72a(^2), 72 J'^2) (straud 2) 



(110) 
(111) 



For simplicity the two strands have the same length Li = L2 = L, so that tj G [0,L], 
^2 G [0, L]. The action is taken to be the sum of the action for strands one and two, plus an 
interaction term between the two strands 



72„72^ 



-3Y.{ll *«/5(^l'^l) + // *a/3(«2,t;2) + // ^i'|(ni,t;2; 
a<j3 V"' "'"l<'"i J J U2<V2 J Ju\,V2 



(112) 



The first line in (112) is the free action for the two strands. In the second line ^^^^(^1,^1) 
and \I/^^('U2, f 2) are the overlap operators (94]) for the two strands 1 and 2. The new operator 
^^^'1(^1, ^2) is the overlap for the contact between strands 1 and 2, defined as 



*S(«i.^i) = *«'(«!, ^2) $i''(Mi,t;2) 



$i'^( 



Ml,f2 



a,b 



(113) 



(114) 
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strand 1 



strand 2 



strand 1 



FIG. 15: The new contact vertex ^]3^{u-:V) and the interaction vertex ^'^'^(n, f) between two 



different strands (labeled 1 and 2). 



and depicted in figure 15 The model can be generalized by assigning different coupling 
constants gi, §2 and gi^2 to the operators '^^^p ^^"^ ■f'o^ simplicity we keep gi = 

92 = 91,2 = 9- 

The two open-strand partition function is defined as 



i 11 E riQ(0)+q2 E ri„(L)+q3 J] r2a(0)+q4 E '^2aiL) 



(115) 



It is calculated in perturbation theory as a power series in g. The term of order g^ is a sum 
of expectation values of products of \E'^^ and operators. Each term is represented 
diagrammatically as a system of two bundles (one for each strand), with a planar system of 
double arches on each strand (for the and \1/^) and of double ribbons between the two 
strands (for the \E'^'^). An example of such a diagram (with only \I'^'^ operators) is depicted 
on figure [161 



q4 



qi 



q2 



12 12 12 12 12 12 

FIG. 16: An example of a two-strand diagram, representing a ^^^i^^'q^^I^^^'q^iIs^js term. Note 
that on the lower strand is left and L right, whereas on the upper strand L is left and right. 



When computing Z^"^^ there is a subtle technical point when dealing with translational 
invariance to factor out the (27r)'^5'^(Sq) terms for each replica. For each replica a one has 
to compute a term of the form 

/ (3)1 . . . <I)2 . . . <I)2 (i)l,2 . . . (f,l,2 i(qiriQ{0)+q2ria(L)+q3r2Q(0)+q4r2c,(-L))\ 
\ ^ /O 
v V V ' 
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Either there is no \I/i_2 operator (-ft'1,2 = 0) and the two strands are decoupled, so that from 
translational invariance we factor out a term 

i2nY5\q^ + qs) (27r)'^5'^(q3 + , (116) 

or there is at least one \l'i,2 operator > 0) and the two strands are coupled, so that 

from translational invariance we factor out a term 

(27r)'^5'^(qi + q2 + q3 + q4) • (117) 

We already note that this subtlety will become fully manifest when under renormalization 
the field r goes to r-y/Z, and correspondingly the momenta q go to q/-\/Z, thus inducing 
different additional powers of Z in (116) and (117). 

We must treat separately the contributions to Z^"^^ according to the number Q of replica 
such that the strands 1 and 2 are coupled. More precisely, a term of order involves Ki 
K2 and Ki^2 with Ki + K2 + -ft'1,2 = K, and may have 2 <Q < 2Ki^2 replicas 
with strands 1 and 2 coupled. Thus we can decompose uniquely Z^'^'^ into the contribution 
of these "Q-sectors" , Z^'^''^\ as 

2^'Hqi,q2,q3,q4;^) = [(2^)''^'(qi + q2)'^'(q3 + q4)]""'^ 

Q 

X [(27r)'^5'^(qi + q2 + q3 + q4)]'^i(''^Hqi,q2,q3;^) • (118) 
The Q = term is the disconnected contribution 

Z^^^'\cir,ci2,ci,-.L) = Z^^\ci,-L)Z^^\ci,-L) . (119) 
Since the interaction is a 2-replica interaction, the Q = 1 term is zero 

^^'''Hqi,q2,q3;^) = o. (120) 

The details of the calculations of these two-strand partition functions will be given in section 
|Vj To identify all necessary renormalisations, we only need to compute the Q = 2 sector 
contribution 2^^'^^ at zero external momenta q = 0. 



IV. UV DIVERGENCES AND RENORMALISATION 

A. Introduction: dimensional analysis 

The LW field theory suffers from short-distance (UV) divergences when po, taken as an 
analytic regularisation parameter, is greater than or equal to 1. Using 

e = 2po - 2 (121) 

as a Wilson-Fisher expansion parameter, it was shown in [28] that at one loop (first non- 
trivial order in g), these UV divergences appear as poles in 1/e, and can be absorbed into a 
renormalization of the coupling constant g (strength of the disorder) and of the strand length 
L. In this renormalization framework, for e > 0, the one-loop calculation shows that there 
is a physically relevant non-trivial UV fixed point g* = 0{e) (with g* > 0), which controls 
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the continuum limit of the LW theory. This UV fixed point separates a weak-disorder phase 
{g < g*), where disorder is irrelevant at large scales, from a strong disorder phase {g > g*), 
where disorder is strongly relevant at large distances. 

In this section, we consider the UV divergences in the RW representation of the model. 
We show that the divergences can be analyzed via the short-distance behavior of the RW 
model, through a multilocal operator product expansion (MOPE). This MOPE is similar to 
the MOPE for polymers (the Edwards model) and for self-avoiding polymerized membranes 
(SAM), and is a generalization of the well-known operator product expansion (OPE) for 
local quantum field theories. 

The RW model in d dimensions is equivalent (for closed RW) to the LW model with 
analytic regularisation where the contact exponent is po = d/2. Hence we denote 

e = d-2. (122) 



Using units of the strand length L, in the action S given by (91-96) the bare scaling di- 
mensions of the base position t, of the position vector r and of the auxiliary fields 7 and 7 
are 

[t] = l , [r] = l/2 , [7] = [7]=0. (123) 

The dimensions of the elasticity r^, of the contact operator $, of the overlap operator 
and of the coupling constant g are 

[P] = -l , [<!>] = -d/2 , [^ = -d , [g]=d-2 = e. (124) 

If the theory behaves as an ordinary quantum field theory, it is natural to expect that 
it has UV divergences for d = 2, is perturbatively renormalizable for d = 2 and non- 
renormalizable for d > 2. This will be true if the short-distance singularities are proportional 
to the operators already present in the theory, and if no new terms are generated under 
renormalization. 



B. The MOPE 

The short-distance singularities can be analyzed via a multi-local operator-product ex- 
pansion (MOPE). The importance of this MOPE was already recognized by Lassig-Wiese, 
(see especially HH]). 

The fact that the short-distance singularities for products of \1/ operators in our RW 
model are described by a MOPE is easy to understand, without much explicit calculations. 
Indeed, the operator 'ifai3{u,v) is a product of two bilocal contact operators $Q,(M,f) for two 
independent replicas a and P 

\1/q,/3(m, f ) = $ Q,(ti, f)$^(tt, f) . (125) 

Each contact operator is the product of the standard bilocal contact operator S(-u, v) for the 
plain RW, times two local operators T(u) and T{v) for the auxiliary fields 7 and 7 

a,b 

with 

E{u, v) = ~6\r{u) - r{v)) = {Anf' 6'{r{u) - r{v)) , T,,(u) = 7a(w)76(w) • (127) 
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It is thus sufficient to analyze the short-distance behavior separately for each replica a; and 
further separately for the $q,, functions of r^, and for the local operators involving only 
the auxiliary fields 7^ and ja- 

When dealing with open strands, one must be careful to take into account the boundary 
operators at m = and u = L. One has to write the MOPE at the boundaries for products 
involving a boundary operator and bulk operators, when some of the points go to the 
boundary. As we shall see, boundary operators are important for the renormalization of the 
model. 



1. MOPE for the r operators 

The MOPE for the S operators is nothing but the standard MOPE for the operators 
of the Edwards model for a SAW, and of the general D > 1 polymerized (or tethered) 
membrane (SAM). This MOPE was studied extensively in[29llMl|3ll[32l[33lElll35l|39]. It 
is obtained by expressing the contact operator as 



E{u,v) = I ^ e''i['-(")-'-(")l . (128) 



TT 



d/2 



One then writes the short-distance expansion of products of vertex operators exp(iqr(M)) in 
terms of normal products, and the short-distance OPE for the massless free field in ID (the 
quantum free particle), using the explicit form of the ID propagator. 

The ID massless propagator Go{u,v) ~ {r{u)r{v))o for the scalar field is the solution of 

= ^^^^^ 

with periodic b.c. for closed strands, and Neumann b.c. for open strands. The 1/L term 
takes care of the zero mode. The propagator is explicitly, for the closed strand 

G^^"''''^{u, v) = {\u-v\- L/2f/L (130) 

and for the open strand 

G°P^°(u, v) = {{u - L/2f + {v- L/2f - L\u - v\) /L . (131) 



We only need the difference propagators, see equation (73), 

{\u — v\ for open b.c. 
,/r , (132) 
\u - v\{L - \u - v\) , , 
for closed b.c. 
Ij 

The product of multilocal operators has a short-distance expansion in terms of other multi- 
local operators, of the general form 

OM.nl-)0,{ulnl...)---0,{ul,ul...) = }) 0(K}) 

(133) 
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This expansion generates an algebra containing all multilocal operators of the form 

0(si, ■ ■ ■ , Sfc) = y d% n (^""^^(50) V'^^5^(r(s,) - ro) j , (134) 

where d and V are the derivative in internal and external space. 

d d 

The coefficients C^^' sf}) of the MOPE are homogeneous functions (or rather 

distributions) of the relative distances -u*— sf . Except when sf is a boundary point, they do 
not depend on sf , but on the differences of the m*. 

For k = 1 (local operators) the most relevant operators (with the highest canonical 
dimension) are 

1 and drdr = (r)^ (136) 

with canonical dimension and —1 respectively (remember that the dimensions of s and 
r are 1 and 1/2). For k = 2 (bilocal operator) the most relevant operator is the contact 
operator 

E{s,, S2) = 5"(r(si) - r(s2)) = J d^v^ ¥{v{s,) - fq) ~5\v{s2) - fq) (137) 

with canonical dimension —d/2. For A; = 3 (tri- local operator) the leading operator is 

H(^)(^i, ^2, s,) = ~5\v{s,) - v{s2)) 6\v{sr) - v{s,)) (138) 
= j d^o 5'(r(si) - ro) 5''(r(s2) - ro) 6''{r{s,) - ro) (139) 

with canonical dimension — etc. 

We give here the explicit form for the MOPE of the relevant operators at leading order, 
as calculated for instance in [321 |33l [39] 

S(mi,M2) = |Mi-M2r'/'l(w) - 7l«i-W2r""/'r(«)' + ••• (140) 
with M = (mi + ^2)72. 

S(Mi,i;i)E:(u2,t;2) =„ (|mi - M2I + l^^i - ^'21)^"*^^ 2(u, t;) + ■•• (141) 



^1 — 



with «=(«! + ^2)72, = (t>i + f2)/2. 

H(-u, t;) S(mi,M2) = |mi - ■U2r'^^^ S(m,'u) + ■•• (142) 



r\ui)r\u2) = 8 6{ui - U2) r\u) + ■ ■ ■ (143) 

with u = {ui + U2)/2. 

Eim,U2) i'iu) \u,-U2\-''^' i\u) (144) 



uo — '■u 
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And for the boundary operator 

S(mi,M2) lb(0) =, Imi-msP'/' lb(0) + ••• (145) 



There is another potential term, but it vanishes because of the Neumann boundary condition 
for the open RW 

r\u) lb(0) = + ■■• (146) 

Similar MOPEs can be written for the product of three or more operators. The coeffi- 
cients C of the MOPE are homogeneous functions of the relative distance between the 
points involved. They have themselves a singular behavior when some of the points coa- 
lesce. These nested singularities are also given by a MOPE, with coefficients which have 
themselves nested sub-singularities, etc. This corresponds to the standard concept of nested 
sub-divergences (associated to Zimmerman forests) in the field theory. It is this nested 
MOPE structure which ensures the renormalisability of the self-avoiding polymer and mem- 
brane models studied in[2SlEDlEIlElE3llMlE3ESlE3ESl[Mlllffl An example 



is given on figure 17 



FIG. 17: Example of nested MOPE 



2. MOPE for 7 and 7 

The propagator for the auxiliary fields is a step function. The OPE for the T operators 
is very simple and is exactly given by 



Tab(u)Trd(v 



cd\ 



{V - U) 6bc Tadiv) + 9{U - V) 5ad Tcfe(f ) 



(147) 



Similarly, for the boundary operators of the open strand one has 



7a(0)Tfcc(t^) = 5ab1c{0) 'Ja{L)Tbc{v) = 5aclb{L) 



;i48) 



This ensures that we keep track of the topology (planar structure) of the diagrams at short 
distance. It is clear that while (147) and (148) eliminate some of the UV-divergences com- 
pletely, they "go along" for the remaining ones, thus do not complicate the analysis. 



3. MOPE for $ in the planar limit 

Using (|T26|) it is easy to obtain the MOPE for the dressed contact operators $ in each 
replica sector (we omit here the replica index a for simplicity of notation), and to take the 
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planar large- limit (A^ being the number of "color" indices a for the auxiliary fields 7a 
and 7a) to obtain the MOPE for the planar pairing operators. At leading order this MOPE 
involves only $, and the local operators 1 and r^, or rather their "dressed versions" 



r(u) 



" j.2 )5 I 



U) 



j;^^la{u)-ia{u) 1{U) = J^^^aa{u)l 

a a 



u 



u) r (u) 



(we shall omit the dressing " " for local operators in the rest of this section). 
The first MOPE (|T40| is unchanged 



$(M1,M2) 



Ml — M2 



1 



U 



Ml — U2 



The two mope's involving two $'s become 

$(Mi,Mi)$(m2,M2) ^^=_ [e{u2-Ui)0{Vi-V2) +e{Ui-U2)e{v2-Vi)] 



X (|mi - U2\ + |mi - M2I 



-d/2 



$(m,m) + 



and 



<I>(m,m) <I>(mi,M2) = [9{u-U2) + 9{ui-u)] IM1-M2I '^^^ ^{u,v) + ••• 



— >u 



(149) 
(150) 



(151) 



(152) 



(153) 



The remaining MOPE's (143 164) for the bulk and boundary operators are unchanged at 
this order. 

One important remark is in order, about the MOPE's (142) and (153). For the simple RW 
contact operator S, (142) is just the product of the MOPE for E 1, times an independent 
S, which can be viewed as a "spectator". As a consequence, the sum of the diagrams 



u u^u 



V ' 

(154) 

which carry a potential UV divergence in the SAW model, are canceled by the counter-term 
for the leading divergence in (151), and no counterterm is associated to these diagrams 
in the Edwards Model. However the corresponding MOPE for the planar pairing model 
is different, since the term associated to a non planar configuration (the second diagram 
in (154), is absent. As a consequence, and as we will see later, the MOPE (153) gives 
a non-trivial UV singularity in the replica-interacting model, and requires an additional 
renormalization , which is absent in the SAW model. 



4. MOPE for 

It is clear that there is an analogous MOPE for the "^a/s, since these operators are products 
of $ operators on two independent replicas a and (3. More generally, one can write a MOPE 
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for any product of multilocal-multireplica operators of the form , 

0.,.,...«,({.}}) = Olds]}) Oliis'^}) ■ ■ ■ Oliis'^}) 



(155) 



where each OJ^^ is a muhilocal operator O* of the form (134) for a single rephca aj. Of 
course, the MOPE is non-trivial only if some of the replica indices a are common to the 
different multilocal-multireplica operators. We give here the MOPE for the \l/ operators 
which are of interest for the renormalization of the model at one loop. For a single \1/ we 
have 



^af3{ui,U2) = \ui ~ U2\-'^ l{u) - - \ui - Uil^"'^ [rliu) + rl{u)] + 



(156) 



+ 



+ 



+ 



U 



For two \& we have 

^a/3(Ml,fl)*a/3(M2,t^2) „^=_„ [9{U2 - Ui)9{vi - V2) + 9{ui - U2)0{v2 - Vi)] 

X {\ui-U2\ + \Vi~V2\)~'^ '^af3iu,v) + ■■■ 



+ 



(157) 



+ 



"1 "2 



^2 ^1 



and 



'^ap{u,v)'^ap{Ui,U2) ^ = ^[9{U-U2) +9{Ui-u)] - U2\ *a 



un — >-u 



(158) 



+ 



Let us note that the product (157) of two \E' which share only a single replica a generates a 
3-replica operator of the form 



,(3) 



(159) 
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which turns out to be irrelevant. But the product (158) involving 3 replicas is relevant 

'^al3{u,v)'^a'r{Ul,U2) ^ = ^ [9{u - U2) + 9{ui - u)] |Mi - M2 f"' (m, f ) + ■■■ (160) 



u u 



1 2 



For the bulk operators and local boundary operators the MOPE is similar to (143)- (164). 
For the bulk we have 



rl{ui)rl{u2) = 8 6{ui - U2) rl{u) + 



(161) 



^ — X- 



and [39] 



— 



-AdQ{ui <u <U2) \ui - U2I ^ l{u) 
+ + [d- 1)6(mi <u < U2)] \ui - U2I 



+ 



(162) 



[il{u) + vl{u 



u 



+ 



The four terms in the square bracket in (162) are obtained by not contracting ia{u) (first 



term), contracting twice (second term), or once (third and fourth term). 
For the boundary operator we have 



^ap{Ul,U2) lb(0) 



Ui - U2\ lb(0) + 



(163) 



"1 "2 
The trivial one is 



il{u) lb(0) = + 



(164) 



u 



+ 



This MOPE structure for the multilocal-multireplica operators can be shown to hold at 
higher order, and to have a nested structure for its sub-singularities, as for the single replica 
case. It is this nested structure of the MOPE which ensures that the model is renormalizable. 
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C. Renormalizability 



We can now analyze the short-distance singularities of the full interacting theory with n 
replicas for e close to 0. The analysis is very similar to the general proof of the renormalis- 
ability of the SAW and SAM models [23 EOl ED E21 E3] , thus we can constrain ourselves to 
an outline. 

We compute the partition function and the correlation functions of the model for open 



strands and generic n, as defined in section III, with c? or e as an analytic regularization 
parameter. Using the standard rules of analytic (or dimensional) regularisation, the integrals 
are calculated for (i 7^ 2, treating short-distance divergences via a finite-part prescription. 
Within this framework, UV divergences appear as poles in the complex d plane. There are 
no long-distance IR divergences since the strand length L is kept finite and acts as an IR 
regulator. 

The UV poles at e = are associated to the terms of the MOPE with the correct power 
counting. They give the superficial UV divergences of the theory. These divergences can 
be subtracted by adding counterterms to the action if they are proportional to the original 
operators in the action and if no new terms are generated under renormalization . 

According to the MOPE derived above, the dangerous UV terms come from the MOPE 
for products of operators of the following form: 

Firstly 



l{u)+ K(n) + rJ(n)) 



(165) 



which describes the MOPE for a product of local and bilocal operators into a local operator. 
The coefficient for the relevant identity operator l{u) gives dimension-full UV divergences, 
i.e. poles in the complex e plane for negative values of e = —2/p, p integer. These do not 
give a pole at e = and thus do not require a renormalization in the minimal subtraction 
scheme. Physically, this term gives a dimension-full UV divergence of the form A~^L, where 
A is a physical UV regulator (with dimension of a mass) and L is the strand length. Since 
the strand length L is kept fixed in our scheme, this divergence is the same for all partition 
functions and does therefore cancel in correlation functions. The coefficient for the marginal 



operators 



gives a pole at £ = and will be subtracted by a wave-function 



counterterm proportional to the free action Sq. 
Secondly 



/ // /// 
/ // /// 



(166) 

which describes the MOPE for a product of local and bilocal operators into a single bilocal 
operator. This gives also poles at e = which are subtracted by a coupling-constant 
counterterm proportional to the interaction term iSint in the action. 

Thirdly, for the open strand, there is a divergence coming from the MOPE for the bound- 
ary operator 1^, here located at m = 



r?«) 



1^(0) 



(167) 
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This divergence is subtracted by a new counterterm proportional to the boundary operator 
lb. It is of course not present for closed strands. 

There are no UV divergences associated to the auxiliary fields 'J^'a'i^'a- This is not sur- 
prising since these auxiliary fields are introduced to organize the topological expansion of 
the perturbative expansion and to construct the planar limit for N ^ oo. 

The analysis of the subdivergences and the fact that the UV poles can be recursively 
subtracted in perturbation theory is similar to the one of [291 EDI EH |32l |33l |35l |39] and 
shall not be repeated here. 

Finally let us stress that the arguments for perturbative renormalizability are valid for any 
value of n (the number of replicas) and of (the parameter for the topological expansion 
of perturbation theory). 

We can now study the renormalized theory, and write the corresponding renormalization 
group equations. 

D. Bare and renormalized action and observables 

Renormalizability of the theory means that, when expressed in terms of the renormalized 
field and of the dimensionless renormalized coupling constant g^i, the theory defined with 
the renormalized action S^[rj^; g^, fi] for the open strand with 

Jo<s<L \^ / ^^i^J J0<u<v<L ^ 

with 

*f^(«, v) = {4ny/'6'{rUu) - rUu))S\r^^{u) - r^,p{u)) (169) 

is UV finite when e = 0. Z and are respectively the wave function counterterm and the 
coupling constant counterterms. They will correspond to a renormalization of the operators 
and \E'. fi is the renormalization mass scale, it has dimension L~^. Zi is a new boundary 
counterterm, it is associated to a renormalization of the boundary identity operator 1^. The 
counterterms are defined as a perturbative series in g^^, of the form 

Z = l + J29nMe) (170) 

fc>0 

Zg = l + J29nBk{^) (171) 

k>0 

Zi = J29nCkie) . (172) 

fc>0 

The counterterm coefficients A^ie), Bk{e) and Ck{e) contain poles at £ = which cancel 
the UV poles of the bare theory. For the general case, these coefficients depend on n and 
A^. No renormalization is required for the auxiliary fields. 

As usual, defining the bare field and the (dimensionfull) bare coupling constant gs as 



rB = raZi/2 , g^=g^ZgZ''fi~' , 6 = d-2, 



(173) 
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we can rewrite the renormalized theory as a bare theory, defined by the bare action iSefrB; Qb] 
given by 

SbK, 9.] = J2f (W + - 5^ // ^U^, v) . (174) 

„ Jo<s<L \^ / ^^pJ Jo<u<v<L 

Similarly, for an open strand, this amounts to renormalizing the boundary operator 1^ 
through 

lf = e-^^lf (175) 

When computing the counterterms and UV-finite observables for the renormalized theory 
from the UV divergent bare theory, one must be very careful with the field renormalization 
and the treatment of the translational zero modes. To see this, let us first consider the 
partition function for the closed strand, Z'^^°^^'^{g, L), as defined by (99). Since the (infinite) 
contribution of the translational zero modes 

f Hd'r^ = [{2nr6\q^ = O)]"" (176) 
has already been factored out, this partition function has dimension 

^closed^^.^) _ ^^ydn _ ^^^^^ 



Taking into account the renormalization of the field r given by (173), the relation between 
the bare partition function, computed with the bare action 5b [re; 9b], and the renormalized 
partition function is 

Z'^^'^^'^igB, L) = Z^^^^^'^Qb. L) (178) 
with Z the wave-function renormalization factor. 



The partition function for one open strand Z^^\c{]g,L) is defined by (104), with q the 
external momentum. It has scaling dimension zero 

i<^'(q;(7,L)~l. (179) 

Since the external momentum q is conjugate to the field r, it is renormalized as 

qB = qRZ-i/2_ (180) 

The relation between the bare partition function for a single open strand (computed with 
the bare action SbIyb] Qb], and expressed as a function of the bare coupling Qb and the 
bare momentum qg) and the UV finite renormalized partition function (computed with the 
renormalized action 5R[rR;(7R] and expressed as a function of the renormalized coupling 
and the renormalized momentum q^) is 

Z^'\cIb; Ob, L) e-'-^^ = Z^\qB;9B,L) . (181) 

Note the additional multiplicative boundary renormalization factor exp(— 2nZi) which is of 
course not present for closed strands, and which is important for the calculations. 

Finally we consider the partition function for two open strands in interaction 
2^^^(qi • ■ ■ q4; 5^, L), defined by (115). We have seen that it must be decomposed into 
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terms associated to the contribution of the so-called "Q-sectors" where Q r eplicas inter- 
act amongst the n available replicas, 2^*^^''^)(qi ■ ■ ■ qa; (7, L), defined by (118). Each term 
2(2,Q)(^q_^ . . . j^g^g different scaling dimension 

Z^''^\qi---q3;g,L)c^[rf'' . (182) 



Therefore the relation between the bare 2-strand partition function and the renormalized 
one is 

^'■''^(qB ■■■■,9.,L) Z-^ e-^«^^ = zf^Hqa ■■■;9.,L). (183) 

Note that the boundary renormalization factor is now exp(— 4nZi) since we have two strands, 
hence four end points. 

The argument is a bit sketchy, but can be made more rigorous by taking into account 
the functional measure T>[r] in the definition of the RW model, and using the fact that this 
functional measure is dimensionless for closed strands, but dimensionfull for open strands, 
with dimension 

P^^] closed ^ ^ ^ p[r](°P"°) ~ [r]"^ . (184) 

The difference can be viewed as the insertion of a contact operator 6'^{r{0) — r(L)) for each 
replica in the open ensemble, which closes the strand and thus leads to the closed ensemble. 



E. Beta function and anomalous dimensions 



We use the minimal subtraction scheme (MS scheme), i.e. define the counterterms Z, Zp 
and Zi in such a way that they contain only poles at e = 0, but no finite part with analytic 
terms in e. 



= 1 +^R— ^ + - 



f Oi9l: 
0{9l) . 



(185) 
(186) 
(187) 



The j3g function for the coupling is defined as the variation of the renormalized coupling 



with the renormalization scale /i. Using (173) we get 

dg 



Pgi.9n) 



—e 



9B 



l + ^^ + (2 + .) 



c/logZ 



n -1 



;i88) 



and our short 



The function jSg is the Wilson flow function considered in Lassig-Wiese 
paper [H]. Its derivative P'g gives the scaling dimension of g^^ (in the sense of Wilson, hence 
in units of mass fi). 

Ag = ^Jgig^) . (189) 

Since the renormalized coupling g^^ is the scaling field associated to the bifocal overlap 
operator ^R(M,t>), the scaling dimension of the operator \& (in units of /i) is 

d 



A* 



(igi 



-Pg{gB 



(190) 
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In [H] we considered the dimension of r (in units of length L ~ l//i), XfIs'r), defined as 



d\og'L{g^) 



[1 + i{9r)] 



with 



d\ogZ{g^) 



(191) 



(192) 



Xrign) equals minus the scaling dimension of r considered as a local operator (hence in units 
of mass /i) 

Ar = -Xrign) . (193) 
Finally the scaling dimension for the boundary operator is (still in units of mass n) 



Alb = 71 (^r) = PgigB 



(194) 



These formulas will be used to derive the anomalous dimensions and the renormalization 
-group flow of the RNA model from the two-loop calculation of the counterterms presented 
in the next section. 

It is also possible to study the renormalization of the contact operator ^gj u, v) 
6'^{ra{u) — ra(f )) and to compute its scaling dimension A$. This is done in section 



IV G 



F. MS and MS subtraction schemes 



We use two slightly different subtraction schemes. The first one stated in (185)-(187) is 



the standard minimal subtraction scheme (MS scheme), where the counterterms Z, and 
Zi are chosen such that they contain only poles in e, but no term analytic in e. It is useful 
to make the connection with the renormalization in the Lassig-Wiese formulation [28], see 
subsection lIVHi 

The second one is denoted MS and is defined as follows: Remarking that the relation 
between the bare and the renormalized coupling constant is g^^ = (7i^ZgZ^+^/i~^, its total 
renormalization factor is 

Zg = ZgZ^ + ^ . (195) 

It contains also analytic terms of order e", n > 0. If we choose Zg to have pure poles in e, 
together with Z and Zi, we obtain the MS renormalization scheme where 



,Zl/2 



gs 



gn^gfJ'" 



(196) 



Z and Zi are still of the form (185) and (187) (but with a priori different coefficients 03 and 



and 



Z„ 



bi 2 

1 +^R— 



+ 



+ Oigl) . 



(197) 



The two schemes MS and MS differ by a fin ite re definition of the renormalized coupling 
constant g^^. In the MS scheme the definition (188) for the new beta function (3g reads 



PgigB 



dgn 
dfi 



9B 



1 ^ rflogZ, 
gn dg,. 



-, -1 



(198) 
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G. Renormalisation and anomalous dimension of the operator $ 

We have seen that the local field r and the overlap operators "^ais are renormalized 
and get dimensions Xr and A^. Similarly the bilocal contact operator involving a single 
replica a, must be renormalized and gets an anomalous dimension A$. This dimension plays 
an important role in the analysis of the model, in particular in the "locking" mechanism 
presented in |28] and discussed below. 



1. Contact operator ^ for a single strand 



To compute the anomalous dimension of we must compute correlation functions in- 
volving this operator. We first consider the contact operator between two points u and v 



for a single strand with replica index a. It reads, from (89), 

<l>a{u,v) = 5'^(r„(M) - r^(t;)) . 



(199) 



(For simplicity of notation, we omit the dressing by the auxiliary fields which is necessary 
to make it a planar operator). It is depicted by a single arch as in Fig. 32 Its engineering 
dimension in units of n is 

d . e 



This operator allows to define the "height operator" h{s) for a single closed strand 

h{s) ■= ^{u,v) . 

J J u<s<v 

Since the theory is renormalizable the operator $ can be renormalized as 



(200) 



(201) 



(202) 



where $b is the bare operator (199) involving the bare field Tb, and Z$ is the renormalization 



factor, computed in either the MS or MS scheme. Its coefficients depend on n (the replica 
number) and contains ultraviolet poles in e. The standard RG analysis shows that the 
scaling dimension of $ in the renormalized theory is 



p{.9r) = ^■fign) = A| + 7<i,(5(j, 



7i>(^r) = Pg{9n) 



d log Z$ 



(203) 



where f3g{gB) is the coupling-constant /^-function given by (188). Its value at the UV fixed 
point g^, i.e. at the freezing transition, gives the contact exponent p at the transition 



(204) 



2. Two-strand contact operator 



It is more convenient to consider the contact operator between two different interacting 
(open) strands, in the same spirit as in subsection III F The two strands are labeled 1 and 
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FIG. 18: Contact operator <I>Jj'^ between two strands 



2 and described by the RW and respectively (and the associated auxihary fields). 
The inter-strand contact operator $^'^(m,i;) is 

<!>l;'iu,v) = ~5%r^M-r2aiv)) (205) 

times the appropriate dressing to ensure planarity. This operator is represented diagrammat- 
ically by a single line between the two strands. Repeating the analysis of the UV divergences 
and the renormalization for this operator, it is easy to see that the singularities come from 
the very same terms in the MOPE as for $, that is from the singularities when a bunch of \l/ 
coalesces on the endpoints of This implies that is renormalized exactly in the same 
way as $, namely that the renormalized operator is 

(206) 



1,2 



with Z$ containing the same counterterms as in (202). 



We consider in perturbation theory the following partition function 



?{2) 



2,|3 ^ 



^(5[ri,]+5[r2^]) 




(207) 



ri„(0)=O r2/3{0)=O;/3^ao 

which contains all the diagrams which contribute to the renormalization of $. This partition 
function is the partition function of 2 strands (replica bundles), with interactions between 
replicas inside each strand 1 and 2, no interactions between these two strands (this would 
result in higher replica operators, not $), and constrained to be in contact via the replica 
ctQ. We have omitted the auxiliary fields and the boundary operators Fi and r2 which are 
present in the definition of the partition functions of open strands (103). The constraints 
on the path-integrals riQ(O) = and r 213(0) = take care of the zero modes, except for r2aoj 
whose zero-mode is fixed through ria^iO) = 0, and the operator $^'^(M,t>). 
This partition function is renormalized as 



(2) 



(208) 



(compare with the renormalization of the Z^'^''^^ partition functions in (184), where there 



are Q replica interactions between the two strands). The diagrams and the corresponding 
amplitudes are discussed in section VE[ and the corresponding counterterm Z$ and scaling 
dimension A<i, are calculated in section VE However, we shall now derive directly a scaling 
relation for A<i, at the fixed point g* ^ which implies a scaling relation between the exponents 
p* (the contact exponent) and = C(5'*) (the full dimension of the r field). 
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3. A scaling relation between A$ and Ar 



We remark that can be written as 



Z. 



(2) 



_d_ 



q=0 



(209) 



Indeed each — ^^^^(^^(q) 



q=0 



gives the partition function for a single open-strand bundle 



with a point marked on one of the n replicas. Taking the square gives the partition function 
for two replicas. The factor of takes care of the fact that the replica ao is not summed 

~(2) 

over in the definition of Z ^ . Using the renormalization of the partition function for one 
open strand given in Eqs. nm and flSlf , is UV finite for 

Z$ = Z-^+'^Z^ . (210) 

Z$ is not subtracted minimally, but this is enough to compute the full dimension of $ at 
the fixed point, using (203) 



p(5'r) = A$(5(r) = 1 + 
= 1 + 



1 - ^yfl,(ft)^'°«'^> 



(2-e) xr{gn) 



(211) 



with Xrign) defined in (191). At a fixed point g* we get the scaling relation 

p* + {2-e)x*=2. 



(212) 

Using the scaling relation p* + (* = 2 from Eq. ([s]), this implies the relation between the 
roughness exponent C* of the height h and the dimension of r 



C = (2 - e)x* = (4 - d)x* 



For d = 3 {e = 1) 



X 



(213) 



(214) 



as expected, since h is the height variable, equivalent to the radial coordinate = |r| for 
the random walk. 



H. Relation with the length-renormalization framework of Lassig-Wiese 

In the original one-loop renormalization-group calculation of [28] by Lassig and Wiese 
the random RNA folding model is not formulated in terms of a dressed random walk model 
with auxiliary fields r. The model is defined as a perturbative expansion in the disorder 
strength g in terms of the overlap operator The renormalization s at 1-loop order are 
extracted from the behavior of the free energy at 2-loop order. Since there is no r field in 
this formulation, there is no renormalization factor Z for this field. Also in their scheme two 
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renormalization s are required: the first one is a renormalization of tlie coupling constant g 
as liere, but tlie second one is a renormalization of the length L of the RNA strand. 

In order to compare our results with those of [28j , and to check that our general argument 
for the renormalizability of the model implies the consistency of the renormalization scheme 
of [2H], we must understand the relation between these two renormalization schemes. This 
relation is in fact easy to establish. 

Let us start from our renormalization scheme. The renormalized action is of the form 
(we omit the auxiliary fields, the associated indices a,b, . . ., the replica indices a,(3, . . ., the 
boundary counterterms, and consider only closed strands) 

= ^ / - Qnl^Mf^-' [ [ ^r{u, v) (215) 

^ Jo<u<L J Jo<u<v<L 

This renormalized action can be rewritten as a bare action Sb in terms of the bare field Tb 
and the bare coupling constant g^^ 

= SbIt^] = \ [ hiuf -gB [[ ^b{u,v) (216) 

^ Jo<u<L J Jo<u<v<L 

provided that we renormalize the field r and the coupling constant g as 

rM=^9ny^\iu) , g^ = g^Zg{g^)Z{g^)''fi-' . (217) 
The overlap operator \1/ is defined as usual as 

^afsiu, v) = 6''{rM - ra{v))6''{r0{u) - rp{v)) , (218) 

so that we have 

^B{u,v)=Z{g^)-''^j,{u,v) . (219) 

We now reconstruct the Lassig-Wiese renormalization scheme from our scheme. We can 
consider the strand length L and the base coordinates u G [0,L] as renormalized length 
parameters 

L = L^ , u = u^e[0, Lh] (220) 

and in our scheme the length L is not renormalized, i.e. the length in the bare action Sb is 
the same as the length in the renormalized action Sb.- 

Instead let us construct another renormalization scheme (hereafter indicated by "tilde" 
superscripts) where the coupling constant g and the length L are renormalized, but where 
the field r is not renormalized. We define the base length parameter Mb and the renormalized 
length Lb as _ 

L^ = Z{g^)L^ , u^ = Z{g^)u^ G [0, Lb] (221) 

with Z((7p.) a renormalization length factor (to be determined). Similarly the bare coupling 
constant ^b is defined as 

g^ = g^ Zgig^) (222) 

with Zg{gT^) the coupling constant renormalization factor. The r field is not renormalized, 
the bare field is simply 

TBius) = r^iun) (223) 
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so that 



(224) 



It is a simple exercise to check that if we choose as counterterms 

Zig^) = Zig^)-' , Zgig^)=Zg{g^)Z\g^) (225) 
we have renormahzed the theory in the Lassig-Wiese scheme, since we can rewrite our 



renormahzed theory, given by our renormahzed action in eq. (215), which gives an UV 



finite theory, as a bare LW theory, with the bare action Sb of the form 



5b [fe 



i2 



0<ub<Lb 



b[Ub,Vs) 



(226) 



0<ub<vb<Lb 



We have thus proven that our renormahzation scheme imphes that the Lassig-Wiese renor- 
mahzation scheme exists and is consistent at all orders. The reader interested in the precise 
correspondence between the one-loop calculations of [28] and our formalism may check that 
in [2H] the calculations are done in a minimal subtraction scheme where the renormahzation 
scale is set by the renormahzed strand length. This corresponds in our formalism to set 
fi = and to chose a minimal subtraction scheme where the counterterms for the cou- 
pling constant g, that is Zg, contain pure poles in e. The relation between the /j-f uncti ons 
is obtained as follows. Insert = = Ljj/Z{g^) = L^Z^g^) in (222), and use (225) to 
obtain 

g^ZgZ'-^%g^)Ll . (227) 



9e 



Comparing (217) and (227), we conclude that 



f ^ 



d 



9B dfi 



(228) 



9b 



The l.h.s. is the /^-function denoted (3{u) (with u = g^i) in LW [28], whereas the r.h.s. is 
our /^-function Pg{gB.) in eq. (188). According to (225) this corresponds to our minimal 
subtraction scheme MS. The precise dictionary between our notations and results and those 
of [2H] are given in table [ij 



this work 


^ by LW 


denomination 






renormahzation scale 


9r 


u 


dimensionless renormahzed coupling 


9b 


90 


bare coupling 


Lb 


Lo 


bare length 


L = L^ 


L 


renormalized length 


= Z 


Zl 


length renormalization factor 






coupling renormalization factor 




/3{u) 


beta function 



TABLE I: Correspondence between our renormalization scheme and the LW scheme |28] . 
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V. 2-LOOP DIAGRAMS 

A. Presentation of the calculation, fixed-length scheme 

This section is devoted to the exphcit calculation of the 2-loop diagrams which are re- 
quired to evaluate the renormalization factors and the renormalization-group functions at 
second order. The main calculation will be performed for open strands in the fixed-length 
framework which is presented and discussed in the previous sections. For completeness and 
comparison with the calculations of [28] we shall also discuss the so-called grand-canonical 
scheme where an integration over the length L of the strands is performed. This scheme has 
the advantages that it allows for an independent check, and that some integrals are much 
simpler. 

We first recall the definition of the amplitudes in the fixed-length framework. The parti- 
tion function for p strands, ^(^^(q ■ ■ ■ ; g^, L) is decomposed into a sum of the contribution of 
the Q se ctors where Q is the number of interacting replicas between two strands, according 



to (115). Each sector gives the partition function 2*^^''^)(q ■ ■ ■ ]g^^L) which is a function of 



the 2p—l independent external momenta q. The 2n — Q delta functions for the conservation 



of external momenta are already extracted and taken into account in (115). To extract the 
necessary renormalization factors, in fact we only need two cases: First, p = 1, and Q = 0, 
with the partition function denoted Z^^\c[; g^^, L). And second p = 2,Q = 2. Note that 
p = 2,Q = reduces to the previous one since 2^^'°) = Z^^'^Z^^^ and the case p = 2,Q = 1 
is trivial since Z^'^'^^ = 0. 



1. One-strand diagrams 

The perturbative expansion of the one-strand partition function Z^^\q;, g, L) is written 
as a sum over one-strand planar diagrams Si 

oo 

Z('\q;g,,L) = Y,9' E c{S,) E,,{ci; L) . (229) 

diagrams 

c{Ei) is the combinatorial factor of the planar diagram £i. This combinatorial (or symmetry) 
factor depends on n (the number of replicas). Ei{q;, L) is the amplitude of the one-strand 
diagram Si. It is given by an integral over the position of the internal arch points of ampli- 
tudes of the general form given by (107). For a diagram of order k (with k double arches), 
by homogeneity this amplitude is of the form 

Ei(q; L) = L'=(2-2e)^.(q^i/2) ^-nq^L ^ ^230) 
where for historical reasons and to make an easy comparison with [28] we keep the notation 

»=^ = l + |. (231) 

Each dimensionless amplitude -Ei(q') depends only on the dimensionless momentum q' = 
qL^/^ and of course on 9. The UV divergences appear in our dimensional-regularisation 
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scheme as poles at 6 = 1 (that is e = 0). Power counting shows that the primary divergences 
(needed to compute the RG functions) are contained in the first two coefficients of the 
amphtudes Ei{q) in a small-momentum expansion in q around q = 0. We denote these two 
coefficients Ai and Bi 

£;,(q) = A, + q25, + 0(q^) , (232) 

(they depend on 6). Finally, since we can express everything in terms of dimensionless 
quantities, we set 

L = 1 (233) 

in the rest of the calculations in the fixed-length scheme. In this case, the calculation of the 
amplitudes -E'i(q) and of the coefficients Ai and Bi for each diagram Si involves integrals 
over the relative successive distances between arch end points on the strand of the form 




Setting = 0, and U2k+i = 1 and denoting vj 
ordered end points we denote such integrals by 



/'U2k — l 

/ du2k ■ ■ ■ (234) 
'o 

= Mj+i — Uj the relative distance between 



(235) 




2. Two-strand diagrams 

The perturbative expansion for the two-strand partition function Z^'^''^\q - ■ ■ ,g,L) is 
written as a sum over the connected planar two-strand diagrams Ci with Q = 2 replica 
interactions between the two strands, 

oo 

Z^'''\ci---;g,,L) = Y,9' E c(C.) ^(q ■ ■ ■ ; L) . (236) 

fc=0 

diagrams 

The two-strand diagrams are now denoted by Ci, and c(Ci) is the combinatorial factor for the 
diagram Ci. It depends on the number of replicas n. Cj(q- ■ ■ L) denotes the amplitude of 
the diagram Ci. It depends on the three independent external momenta q = {qi, q2, qs}, on 
the length of the two strands (which are taken to be equal for simplicity Li = L2 = L), and 
on 9 (or d or e). Following the discussion of the previous sections, the amplitudes Cj({q}; L) 
have a scaling form similar to ( 230[ ), 



a({q}; L) = L'W-o) CK{qLV2}) . (237) 

The UV divergences which are needed to compute the RG functions are contained in the 
first coefficient of the 2-strand amplitude at small momenta q. Similarly we can deal with 



dimensionless amplitudes by setting as in (233) L = 1. Thus we denote the dimensionless 



zero-momenta 2-strand amplitude coefficient for the diagram Ci as 

a=a({q};L)|q^o,L=i • (238) 



These coefficients still depend on 9 (or d or e) and have UV poles at 6* = 1. 
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B. Grand-canonical scheme 

We had defined above the canonical, i.e. fixed-length, partition function for open strands, 
(^q . . . -^g^^L^, We now define the partition function for open strands in the grand- 
canonical scheme. This scheme leads to several important simplifications in the calculations, 
and therefore offers a valuable check. 



1. 1-strand diagrams 
We start with the 1-strand diagrams. The partition function is defined as 



POO ^ 

^0 



(239) 



As an example, the free (non-interacting) single-strand partition function for the n-times 
replicated open RNA is 

i(^)(q,L) =e-"'i'^ . (240) 

Thus, its Laplace transform reads 

(q, r) = / dL e-^-Z(^) (q, L) = — ^ . (241) 



According to (230), the contribution of diagram % to the 1-strand partition function at order 
is of the form 



Its Laplace transform reads 



(242) 



dLe-^^E,(q;L) 



1 



r + raq 

Thus we can write for diagram i at order 



[r(l - ke){T + nci^ f'Ai + q2r(2 - fe)(r + nci^ f'-^B, + . . . ] (.243) 



El{q,r)g 



[g^{r + nce)T{l-e)t 



r(l - ke) V{2-ke) B,q' 

Lr(i -£)'=' + 



r(l -e)'^ r + nq2_ 



(244) 



To simplify calculations, in the next subsections we will evaluate all amplitudes with the 
combination t + nq^ = 1 kept fixed. We define and BJ to be the coefficients in the square 
brackets of the above expansion, i.e. 



^ - r(i-£)fc^ ' 



p. ._ r(2 - ke) 



(245) 



These are the coefficients which we calculate in the next subsection. Note that the division by 
l/r(l— e)*^ was introduced to have the simplest possible coefficients, especially no derivatives 
of the F-function. This is achieved by replacing g^ by 5fB/r(l — e) in the bare action (174). 
As we will see from the explicit calculation, this replacement leads to identical RG-functions 
up to 2-loop order in the both schemes. 
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2. 2-strand diagrams 

We now discuss 2-strand diagrams, for which the differences are more important. In 
principle, there are two possible definitions: One could either take 2 strands of the same 
length L, and then perform the Laplace transform, or one can take a length Li for the 
first and a length L2 for the second strand, and then perform a Laplace-transform for each 
of the two strands separately, using the same chemical potential. The first method will 
give essentially the same results as in the fixed-length framework, and will not allow for 
any simplifications. On the other hand, the latter method, which we shall adopt in the 
follow ing, leads to important simplifications, see diagram Cn, eq. ( 321[ ) to ( |343 ), and eq. 



(344). Formally, if the contribution of diagram i to Z*^^'^)(q . . . ; (^b; Li, L2) (with length Li 
for strand 1 and length L2 for strand 2) is for the only needed contribution at q = given 
by Ci{Li, L2)gQ, then we define 

-1 /'OO /"OO 

^^^^^fxT^/o "^^Vo d^2e--(^^+^^)a(Li,L2) . (246) 

These are the coefficients which we will calculate in the next subsection. Of course, we 
will not calculate them first for fixed lengths, and then do the Laplace transform, but take 
advantage of the fact that after Laplace-transform many integrals factorize quite nicely. 
Note that we have again consistently divided each diagram by r(l — e)'^, as was done for 
the single-strand diagrams, due to the replacement of by gB/T{l — e) in the bare action 



(174). The coefficients Cl will be calculated in the next subsection. 

To keep notations compact, we will indicate the corresponding integrals in the fixed 
length scheme and in the grand canonical scheme respectively with superscript L (for the 
length-cutoff) and with superscript r (for the chemical potential). Equivalent integrals are 
denoted by "=" to distinguish them from the same diagram, so e.g. Ci = C2 means that 
diagrams Ci and C2 have the same integrals in both schemes, although their combinatorial 
factors, denoted c(Ci) and 0(6*2) may be different. Writing = C|' denotes that this is 
true only for the integral in the fixed length scheme. 



C. Explicit diagrams: single-strand 

1. Order k = 1: 

There is only one diagram. 

£2= /'J (247) 

r s t 

c(^2) = 1 X n(n - l)/2 , ^2(q) = ^ e^^'^^'^'^ (248) 



r+s+t=i Jo ' ' r(3-2^) (2-2^^)(l-2^^) 

= ^ - 1 + 0{e) (249) 

'~Ls^t ~ r(4-2^) -(3-2^)(2-2^)- . ' + ^2^°^ 
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B^^-- (252) 



2. Order k = 2 

There are 5 inequivalent diagrams. 

^3 = / /-v) /'v> / (253) 
r s t u V 

c{8^) = 1 X n(n - l)/2 , E2{q} = I e2(^+")^9^ s-^^^-^^ (254) 

J r,s,t,u,v 



' r,s,t,u,v 

A\ = f s-^'u-^' = \ j s-^'u-''\l -s-u 

J r,s,t,u,v 



2 

1 



S+U<1 



f dxx^-^\l-xf I 

Jo Js+u=l 



2 

^ i r(2-4g)r(3) r(i-2g)r(i-2g) ^ r(i-2g)^ = l_2,,o_:^,, . 

2 r(5 - 46) r(2 - 4^) r(5 -46) £2 £ + l^"^ 6 ^ ^ 

(255) 

Bi= [ 2(3 + u)s-''u-'' = dx x'-'\l -x)' [ s-'' u-'' 

Jr+s+t+u+v=l Jo Js+u=l 

_ r(3 - 4^)r(3) r(i - 2^)r(i - 26) ^ 2(2 - 46) 

r(6 - 46) r(2 - 46) " (5-4^) ^ 

>15 ^ (^2)' (257) 
B; = BlAl (258) 



2 _ + ( -12+ — ) +0(£) (256) 



r,s,t,u,v 



0(84) ^ 1 X n{n - 1) /2 (260) 

yi^ = / (s + «)-^^r2^ = /" {i-{s + u + t)){s + u)-''H-^' 

Jr+s+t+u+v=l Js+u+t<l 

= / xil-{x + t))x-^H-^' - -y) f x'-'H-'' 

Jx+t<l Jo Jx+t=l 

_ r(3 - 4^)r(2) r(2 - 26)r{l - 26) _ r(2 - 26)r{l - 26) 
~ T{5-46) r(3-4^) ~ r(5 - 46) 
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r(-e)r(-i-£) -1 11 



/ 2{s + t + u){s + u)-'h-'' = 2 dyy'-^'{l - y) f x'-^'r 

Jr+s+t+u+v=l Jo Jx+t=l 

r(4 - 4^)r(2) r(2 - 2^)r(i - 29) _ 2(3 - ao) 



(263) 



r(6 - 49) r(3 - 4^) (5 - 49) 

2 6/ TT^ 

2 + 4e 2,2 
+ 



■A 



2 + - + ( 10 - — ) + (264) 



5I = ;^^ = ^ + --2 + 0(5) (265) 

^5 - / ,;--;> / , E,{q)^ [ e'^^-^^^"^%-'' u-'' = Es{q) (266) 

J r,s,t,u,v 

c{£5) = I X n{n - l){n - 2) , A5 = As , B^ = Bs (267) 



(268) 



c{Eq) = 1 X n(n - l)(n - 2) , £;6(g) = / e('^+='*+'")^«' (s + i + u)-^ {s + li)-^-'^ (269) 

J r,s,t,u,v 



A\= I (s + t + u)-'is + u)-'t-^' 



r+s+t+u+v=l 

/ {1- {s + t + u)){s + t + u)-\s + u)-H'^^ 

J s+t+u<l 



s+t+u<l 
1 



dyy'-'\l-y) / x^'V 

^a;+t=l 

_ r(3-4g)r(2)r(2-g)r(i-2g) _ _ A , 9 . ..7.. 

" r(5-4^) r(3-3^) " 46 + 4 + ^^^^ ^'^^ 

r(-2£-i)r(-e- i)r(i- 1) 39,1, , ^ 



= {2s + St + 2u){s + t + u)~^ {s + u)-H 

J r+s+t+u+v=l 

= fdyy'-'\l-y) [ i2x + 3t)x'-H-'' 

Jo Jx+t=l 



r(4 - 4^)r(2) r(2 - ^)r(i - 29) (2(2 -9) + 3(1 - 2^)) 

r(6 - A9) r(3 - 3^) (3 - 3^) 

(3- l^)(7-8g) 
(5-4^)(3-3e) ' 

1 +1+ O(^) (272) 

2£2 2£ V2 24 / ^ ^ ^ ' 
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^6 = 



2^v^(4e + l)r(-2£) _ 1 

■£(£ + i)r(i-f)r(|-|) "2^ 



— — ( 

2e ^ 24 ^' 



-36 + TT^) + 0(e) 



(273) 



^^7 



0(^7) = 1 X n(n- l)(n-2)(n-3)/8 , 



n 2 



^2 ' 



-By = 2A2B2 



(274) 
(275) 



D. Explict diagrams: 2 strands 

We now consider 2-strand partition functions. For simplicity both strands have the same 
length L, or the same chemical potential r, respectively. Note that in order to transform 
from one ensemble to the other, one would need to know quantities for different lengths, 
or different chemical potentials, which are more difficult to calculate. This can be seen as 
follows: If one inverse-Laplace transforms results for fixed chemical potential, one obtains 
results where the sum of the lengths on both RNA-strands equals L, and not their individual 
lengths. 

We recall that there are subtleties for defining and computing multi-strand correlation 
functions, which are discussed above. We only need the Q — 2 diagrams, i.e. the connected 
2-strand diagrams with inter-strand interactions between Q — 2 replicas. It is also sufficient 
to evaluate the amplitudes at zero external momenta q = 0. At two loops we need the 
diagrams with A; < 3 interactions. 



1. Order k = 1 



17 



c(Ci 



1 X 



n{n — 1) 



(276) 



2. Order k = 2 



Co 



^7 



, C2 



-26 



C{C2) = 4 X 



n{n — 1) 



(277) 



r+s+t+u=l 



(3-2^)(2-2^)(l-2^) 

— A'^ 
2—^2 



,.2-2.,L_ r(3-2g)r(i) ^ 
.^.^1 r(4-2^) 

+ 0{e) 



CI = 



(278) 
(279) 
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17 



, Cs= [ r^' = C2 , c{Cs) = 4 X n(n - l)(n - 2) 

J r,s,t,u 

(280) 



C4 = 



17 



c{Ci) = 2 X 



n(n — l)(n — 2)(n — 3) 



O4 _ i^i 



_ r(4 - 2d)_^, 



r(l-£) ^ £(!+£) £ 



- 2 + 2£ + 0{e^) 



7 




7 












/ ..' / 



, C5 



r,s,t 



I {s + v) 

J u,v,w 



-20 



C{C^) = 1 X 



n{n — 1) 



C^^j [ {s + vy^ ^ r ds f dv{l-s){l-v){s + v) 

Jr+s+t=l Ju+v+w=l Jo Jo 



-26 



8-4-^ + 2^-5 1 1 

h 



(2-^)(l-^)(3-2^)(l-2^)2 e 
1 



- 21n(2) 



(281) 

(282) 

(283) 



(284) 



(285) 
(286) 



3. Order k = 3 



/ 






/ 
















7 7 / / 7 



C(C8) = 1 X 



n{n — 1) 



, C8= / / {s-rWr''\t-rX)-''' 

' r,s,t,u J v,w,x,y 



' r+s+t+u=l J v+w+x+y=l 



= / / 

Jr+s+t+u=l Jv-\ 

= / / {l-{s + t)){l-{w + x)){s + w)-^\t + x)- 

J s+t<l Jw+x<l 



26 



(287) 
(288) 
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K{s,x) 



A{i-ef{i-2ef JO 



1 .1 

ds I dx K{s, x)K{x, s) 
Jo 



s'-''' - (1 + s - xf-""' + 2(1 -x){l- e)s 
l^-l + MlG)^ 



1-26 



r\2 



CI ^ {CD 



TT 



2- :_ -41n(2)2 + 81n(2) 



(289) 
(290) 
(291) 
(292) 



17 



: C,= I / is + v)-''x 

' r,s,t J u,v,w,x,y 



26-20 



c{Cq) = 4 X 



n{n — 1) 



L aL 



9 — -^9 ^2 



J r+s+t J u+v+w Jo Jv+w<l 

1 



(1 - 2e){2 - 2e){3 - 26) 



f dv [{1 + vf-^^ - v^-^\2 -29 + v)] (1 - 
Jo 



(293) 

(294) 
(295) 

^3-26* 

(296) 



The last integral is a hypergeometric function, finite when 6-^1. Therefore 
1 



^9 



(l-2^)(2-2^)(3-2^) 
-^-^ + ln(4) + 0(£) 



-(7-4e)r(3-2^)r(4-2^) ^ 2i^i(l,2^-2,5-2^,-l)" 



r(7 - 46) 



(4 - 26) 



= -1 + Q - 2 ln(2)) i + (^-2 + ^ + 2 ln(2)2 - 2 ln(2)) + 0{e) 



<^9 — <^5^2 



(297) 

(298) 
(299) 



•-10 



C 



10 




-20„,-2e 



c(Cio) = 2 X 



n{n — 1) 



: C,o^ I / {S + V + X)-''W 

' r,s,t J u,v,w,x,y 



2e„.-2e 



(300) 
(301) 



[s + v + x)~^"w 



r+s+t=l J u+v+w+x+y=l 



f ds{l-s) [ z{l - z - w){s + z)-^^w-^'^ 

Jo J z+w<l 

^ (1 - Z)-' ((£ - Z)Z-' + {Z + 1)-'Z) 



£2(£ + l)2 



58 




59 



•-15 



17 



c(Ci5) = 2 X 



n{n — 1) 



Ui5 = Ui4 = U2 



(316) 



C 



16 



•-17 



/ 




/ 


I-- 






1.- / 

n(n - l)(n - 2) , 









16 



/ / ' 

J r,s,t J u,v,w,x,y 



28-26 



(317) 
(318) 



Ci7= / / {s + v + w + x)~'^{s + v + x)~'^w 

' r,s,t J u,v,w,x,y 



o„..-2e 



c{Cn) = 2 X n(n- l)(n-2) 



(319) 
(320) 



The integral C^j is hard to calculate. In order to simplify matters, we calculate instead its 
derivative C(y w.r.t. L. Using that Cfy ~ L^~^^, this gives 



C 



n— I ds {L — s) z {L — z — w){s + z + w) {s + z) w 

J z+w<L 

' c: 



e^..-28 



L=l 



r-L 



(321) 



C'lT^ f ds f z{2L-s-z-w)(s + z + w)~^(s + z)~^w 

Jo Jz+w<L 



e„..-2e 



L=l 



5-40 



(E + F + G) , 



(322) 



where we have again taken the derivative w.r.t. L. The different contributions, named E, 
F, and G are 



/ 



E= I z{L- z-w){L + z + w)-\L + z)-'^w 

'o<z+w<L 



o„.-2e 



El + E2 



(323) 



L=l 



which we integrate by part w.r.t. w. The boundary term is Ei, the remaining one E2: 



El 



1 



1-2^ 



dz [z{l -z-w){l + z + w)-\l + z)-^w^-''^] lzl~' = 



(324) 



1 r'^'^ 

E2 = / dz dwz[{l + z + w) + e{l- z- w)] (1 + z)-\l + Z + w)-^-% 

1 - Jo Jo 



i-e„„i-20 



E' 



1-29 



1 t-l-z 

dz / dwz{{l + z) + e{l- z)){l + z)-^-'''w 
Jo 



l-2e„,,l-26» 



(325) 
(326) 
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and E2 = E2 — E2. The decomposition is such that E2 can be integrated analytically over 
w, whereas E2 is finite, s.t. it can be expanded in e. 



E' 



1 



1 - 292 - 29 



dzz{{l + z)+ 9{1 - z)){l + z)-^-''{l - z) 



\2-2e 



4e 



+ 



4 



ln(2) 



+ 0{e) 



- + 0{e) 
The next contribution is 



z{L - s){L + sy%s + z)-% 



0<s<L J z+w=L 



F1 + F2. 



(327) 
(328) 

(329) 



L=l 



Again we split into a contribution Fi, containing the divergent contribution, but integrable 
analytically, and F2, the finite rest to be evaluated at e = 0: 



El 



w 



-26 



[{1- s)w{-s + 9 -l){s + r 



-29-1 



0<s<l ^0<to<l 

4-e-i [_2 {4e^Q - 1) + 1) 9{29 - 3) - 4^ + 2] 
(0(2^ -3) + 1)2 

5 - lnf2l 1 

^ ^ ^ ' [-9 + 21n2(2) + ln(512)] +0(£) 



4 



's - l)s 



0<s<l io<«;<l + -W + l) 



0{e) 



- — [27r2 - 3 (5 - 4 ln2(2) + ln(64))] + 0{e) 



24 

The last contribution is 



G 



-d„..-2e 



2z{s + z + w) %s + zy^w 



0<s<L J 0<z+w<L 



(330) 



(331) 



(332) 



L=l 



which is the most divergent one, since both the integrals over w as over the global scale will 
give a pole in 1/e. To isolate the global pole, we again derive w.r.t. L, reducing this time 
the number of integrations by 1. Since G ~ we obtain 



G 
H 

K 



1 



4-40 



{H + K) 



-d„..-2e 



2z{L + z + wy\L + z^w 



0<z+w<L 



2z{s + L)-\s + z)-^w 



0<s<L J z-\-w=L 



L=l 



L=l 



Hi + H2 + H3 



K1 + K2. 



(333) 
(334) 

(335) 



From both H and K, we need also the finite term, since it is multiplied by 1/e from (333) 



but we will even calculate the next sub-leading term. We use again the technique to split 
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H into parts which contain the divergent contributions Hi and H2 and which are doable 
analytically, and a part for the rest, which is finite and can be expanded in e. 



Hi 



l-z 



20„,-20 



dw 2z{L + z)~^'^w 



Ho = -29 / d^ 



dz dww^-^'^ziz + iy 
Jo Jo 



20-1 



' 2{z - l)z {1 - z') 
29-1 

e + 2 '■^ 



2\-2e 2 - 



0rr(-6) 



2(s + l) 



f d^(l-^)-=^(^ + l) 
^0 



(336) 
-'-^ (337) 



Since the last integral is convergent, we can expand it in s 



Ho 



e + 2 
e 

e + 2 
e 

1 



dz 



(^ + 1)3 {z + iy 

1 3e 7e^ o 
- + — + — + Oie^ 
8 16 32 ^ 



ln(i=f)£ ^ln2(l=f)£2 



2(^ + 1)3 



H. 



dz 



l-z 



dww~'^^z{z + 1) 



-20 



{z + l)\w + z + l)-' + 



w9 

V+1 



- 1 



(338) 
(339) 



Having calculated separately the most diverging terms in w, Hi and H2, the remainder H3 
is finite, since the square bracket is of order w'^. Thus the whole expression can be expanded 
in e, and then integrated analytically. This is straightforward, but since the intermediate 
expressions are rather cumbersome, we only give the final result: 



H. 



8 



+ ln(2) 



+ 



llTT^ 

96 



-(9-ln(2)(18 + ln(2))) 



e + 0{e') 



(340) 



We now treat K, which we split into a diverging integral Ki, doable analytically, and a 
convergent rest K2, which we expand in e: 



Ki 



w-^' (2(s + l)-2«- 2(5 + 1)- 

0<s<l JO<w<l 

]^ _ n _ 2-£-2^ I o^ Q _ n-£-2 



-20-1 



w 



's-9 + 1)) 



e{e + 2)-3- 2" 



e(e + l) 



= — + 



Ko^ 



4e 
2 



31n(2) 



+ 



3-21n(2) 



31n2(2) 



20„,-20 



0<s<l JO<w<l 

1 ln(2) 
8 T~ 



+ 



s + l^w 
■^_23.^ 



{s + ifiw -l){s-w + 1)-^ + w 



w9 

V+1 



96 



+ ln(2) - 



71n2(2) 



e + 0{e^) , 



(341) 
1 

(342) 



where in order to integrate K2, we have again expanded in e, which is justified since the 
square bracket is of order w^. 

Taking all terms together, the final result for Cf^ is 



5 



ln(2) 



+ 



4 + ^-lln(2) + ln(2)^' 



+ 0{e) 



(343) 
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The diagram in the grand-canonical scheme is much simpler, and all integrals can be done 
analytically, using the standard Schwinger representation. The result is: 



2-(^+i)V^csc (f ) r - e) _ _^ ^ ^ ^ _ ^ ^^^^ 



£2(£ + l)r(l- f)r(|-l) 4£2 ' 8£ ' V 8 48 

The remaining diagrams are: 



(344) 



Cl8 — 



c{Cis) = 4 X n(n - 1) (n - 2) , Cis = C 



•-19 



:7 



C20 — 



-21 



•-22 



c(C2i) = 4 X n(n- l)(n-2) 



5-4^ 



3 3 
'4£ 4 



+ 0{e) , = Al 



13 



(345) 



c(Ci9) = 4 X n(n - 1) (n - 2) , C19 = Cn 



(346) 



c(C2o)=2xn(n-l)(n-2) , C20 = C; 



12 



(347) 
(348) 
(349) 



c(C22) = 4 X n(n - 1) (n - 2) , C22 = Q 



21 



(350) 



^23 — 



-24 



<^24 



c(C23)=4xn(n-l)(n-2) , C23 = Q 



c(C24) = 4 X n(n- l)(n-2) 



J r.s 



21 



(351) 
(352) 
(353) 
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C25 — 



c(C25) = 4 X n(n - 1) (n - 2) , C25 = C, 



24 



17 



(354) 



26 



c(C26) = 4 X n(n - 1) (n - 2) , C26 = C; 



24 



-ill tL_ 



Zl 



17 



(355) 



C27 — 



c(C27) = 4 X n(n - 1) (n - 2) , C27 = C; 



12 



:7 



(356) 



-28 



c(C28) = 2 X n(n - 1) (n - 2) , C28 = C; 



12 



:7 



(357) 



C29 — 



c(C29) =4 X n(n- l)(n-2) , C29 = C; 



14 



\ — ~ _ I. ■ 



:7 



(358) 



-so 



c(C3o) = 4 X n(n - 1) (n - 2) , C30 = C^a 



(359) 



C31 = 



c(C3i) = 2 X n(n - 1) (n - 2)\ C31 = C; 



14 



(360) 



C32 — 



c(C32) =2 xn(n-l)(n-2)2 , C32 = C; 



14 



(361) 



33 



c(C33) = 2 X n(n - l)(n - 2f , C33 = C 



Cm — 



-35 



c(C35) = 2 X n(n - l)(n - 2)(n - 3)/4 , C^, = C^A"^ 



t-35 - p^^ 



2-20 
2 



(£ - 2) (2r(3 - ef + (e - 2)r(3 - 2e)) 
2(e2 + e)2r(3-£)2 



14 



(362) 



c(C34) = 2xn(n-l)(n-2)2 , C34 = C; 



14 



(363) 

(364) 
(365) 



(366) 



C36 — 



Z7 



C(C36) = 2 X 



n(n — l)(n — 2)(n — 3) 



^ .L ^ g3^6r(2 - 2S) _ l-2£ _ 1 _ 4 



(367) 



-- + 7 + C>(£) (368) 



■^37 



17 



C(C37) = 2 X 



n(n — l)(n — 2)(n — 3) 



.r _ C'3^7r(2 - 2£) _ -1 + 2£ _ 1 3 



Cl, = C{A\ , CI, = ^p(71^)7' - ^2(; +7) " " ^ + ^ ~ ^ + ^^^^ 



C38 — 



C: 



38 



J r,s,t,u,v 



s + t + u)-%s + u)-h-^\r + s) (371) 



(369) 
(370) 
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c(C38) = 4 X n(n - l)(n - 2){n - 3) 

Jr-i 



^38 



' r+s+t+u+v=l 

Symmetrizing {r, s} {v,u} gives: 



(s + 1 + + u)-'r''\\ - t)/2 



r+s+t+u+v=l 
1 



-0-64.-20 



x{i-t- x)(i -t)(x + ty^x'^t 



x+t<l 

1 



1 



1-04.-26 



dyy'-''\l-y) / -{l-yt)x'-'t 

Jt+X=l ^ 

1 r(3 - 4^)r(2) [T{2 - e)T{i - 20) (3 - ie) r(2 - 9)t{2 - 2^) 

(5 - 46) r(4 - 3^) 



2 r(5 - 49) 
_ 1 _ 7 

~4^ ~ 8e 
r(2 - 2e)C. 
r(l -£)2 4^2 8£ ' V8 48 



r(3 - 3^) 



38 



1 3 /15 7r2\ , 



(372) 

(373) 



(374) 
(375) 



C39 — 



C39 = C21 



c(C39) = 4 X n(n - l)(n - 2)(n - 3) 



(376) 



Cac) — 



C40 = C21 



c(C4o) = 4 X n(n - l)(n - 2)(n - 3) 



(377) 



t 

C41 = C12 



:7 



c(C4i) = 2 X n(n - l)(n - 2)(n - 3) 



(378) 



'42 





/ 










/ 



C(C42) = 4 X 



n{n — l)(n — 2)(n — 3) 
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= 



r(2 - 2e) 



(5 -4^)(2- 2^)2(1 -2^)2 £2 
r(2 - 2e) 12 



3 + 0(£) 
(379) 



TT 



+ ( 3 + y ) + (380) 



C43 — 



C(C43) = 4 X 



n(n-l)(n-2)2(n-3) 



_^L/lL — f^T AT 



(381) 



C44 — 



c{CJ = 1 X 



n(n- l)(n- 2)2 (n- 3)^ 
8 



r.r _ . ,L^2 ^(4-2g)2 ^ _ 4 
r(l-£)2 £2 £ 



C44 — (^2 )^ ) C'Xi _ {A2 



+ 8 + 0{e) 



(382) 



Ca5 — 



C(C45) = 4 X 



n{n- l)(n-2)(n-3)(n-4) 



C45 = C^A2 , CJg = 



r _^L.Lr(6-4^^) _ 1 3 



r(l-£)2 £2 £ 



TT 



+ — + 4 + 0{e) (383) 



46 



37 



C(C46) = 2 X 



n(n - l)(n - 2)(n - 3)(n - 4) 



_ .L r(6-4g) ^ 1 _ 4 

^r(l-£)2 £2 £ 



C46 = C^A^ , CJg = /I3 



- - + 7 + 0(£) 



(384) 



'47 



37 



C(C47) = 2 X 



n(n- l)(n-2)(n-3)(n-4) 
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3 15 , 

■- + — + £ 

As A 



Z7 



■^48 



C(C48) = 2 X 



n(n - l)(n - 2)(n - 3)(n - 4)(n - 5) 
16 



(385) 



(386) 



Cap — 



:7 



C(C49) = 4 X 



n(n-l)(n-2)(n-3) 



_ r^.L _ ,L,. r(4-2g) _ r(2-2e) . _ 



1 3 



r(l-£) ~ r(l-£)2 £2 £ 



+ e»(£) (387) 



■'SO 



c(C5o) = 4 X 



n{n — l)(n — 2)(n — 3) 



C'so = (^49 



(388) 



Zl 



17 



C51 — 



c{C5i) = 4 X n(n - l)(n - 2){n - 3) 



(389) 



Ck2 — 



^52 — ^2 



— ^2 



c(C52) = 4 X n(n - l)(n - 2){n - 3) 



^4-4. / ,-2. 3-2. _ ,L 1 r(i-2^)r(4-2^) 

.-^.^Z -^^5-4^ r(5-4^) 
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'-^52 



1 7r2 

E 6 



rfi-£)2 



1 3 



+ 5 + 0{e) 



(390) 
(391) 



53 



c(C53) = 4 X n(n - l)(n - 2)(n - 3) 



,-}~ - — ; L_ 



C53 = C*!! 



(392) 



54 



c(C54) = 2 X n(n - l)(n - 2)(n - 3) 



C54 = C12 



(393) 



E. Diagrams and amplitudes for $ 

We compute in p ertur bation theory up to two loops (order g^^ the partition function 
Z^\g,L) defined in (207). g is the (bare) couphng constant and L the length of the open 



strands. It is written as a sum over the diagrams (labeled by V^) involving two-strand 
bundles with a single contact line (the $ operator) between the two strands and arbitrary 
(planar system of) double arches within each strand bundle (the \Ef interaction operators) 
The perturbative expansion thus reads 



4') (^, L) = 5^ / ciV.) A (394) 



?(2) 

k=0 V, 



T>i are the diagrams, c{T>i) their symmetry factor and Di the amplitude. The amplitudes 
Di can be deduced from the amplitudes Cj already calculated for the 2-strands diagrams Cj 
considered for the renormalization of g. 

At order k = there is only the trivial diagram 



Vi = 



c(Vi) = 1 , Di = l (395) 
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1. order k = 1 



At order 1 there are 2 diagrams 



c{V2) = 4(n - 1) , D2 = Cs (396) 



c(D3) = (n - 1) (n - 2) , ^3 = C4 (397) 



order k = 2 



At order 2 there are 19 diagrams 



ZL 



c{V^) = 4(n - 1) , ^4 = Ci9 (398) 



^4-^1^^ — ; — ^ 



:7 



c(P5) = 4(n - 1) , L'g = C18 (399) 



c{Ve) = 2(n - 1) , Dq = C20 (400) 



c{Vr)=4{n-l){n-2) , D, = C53 (401) 



c(p8) = 4(n - 1) (n - 2) , L>8 = C51 (402) 
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c(P9) =4(n-l)(n-2) , Dg = Q 



52 



(403) 



c(Pio) = 4(n - l)(n - 2) , L>io = C40 (404) 



c(Pn) = 4(n-l)(n-2) , = C39 (405) 



c{V^2) = 4(n - l)(n - 2) , D^2 = C38 (406) 



c{Vis) = 2(n - 1) (n - 2) , ^13 = C54 (407) 



14 



c(Pi4) = (n-l)(n-2) , 



36 



(408) 



c(Pi5) = (n-l)(n-2) , 



15 



37 



(409) 



16 



c(Pi6) = 4(n - 1)2 , 



Die = C31 



(410) 
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V 



17 



V 



18 



2^19 = 



n 



20 



p. 



21 



p. 



22 



, c(Di7)=2(n-l)(n-2)(n-3) , 



, c(Pi8) =2(n-l)(n-2)(n-3) , 



, c(Pi9) =2(n-l)(n-2)(n-3) , 



, c(P2o) =2(n-l)2(n-2) , 



, C{V22) 



{n-lf{n-2f 



D 



17 



45 



(411) 

Dl8 = C46 

(412) 

-D19 = C47 
(413) 

(414) 



^ (n-l)(n-2)(n-3)(n-4) 



(415) 

D22 = C44 
(416) 



VI. RG-FUNCTIONS AND RESULTS 

A. Calculation of the UV poles and finite psirts at 2-loop order 

In the last section, we have computed the first two terms in the small-momentum expan- 
sion (in powers of the external momenta q^) of the bare (unrenormalized) single-open-strand 

partition function Z^j^\cIb; Qj^, L), as well as the two-open-strand partition function in the 

'^(2 2] 

Q = 2 replica sector Z^' [q^si-, g-si L) at zero momenta (q^j = 0). We consider here the 
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fixed length L ensemble and normalize the strand length L to 

L = 1 (417) 

to simplify the calcula tions . This is equivalent to deal with the dimensionless momenta 
q' = qL^/^ (see section VA). To order and order q| this single-strand partition function 
reads _ 

Z^ji\q,;g,)=A{g,) + qln{g,) + 0{q,^) (418) 

with 

A{g^)=Ao + g^Ai + glA2 + 0{gl), M{g^) = Mo + g^M^ + glM2 + Oigl) (419) 
The two-strand partition function is 

Z^B''\^Bi, 9b) = C{g^) + 0(q/) (420) 

with 

C(^?b) = ^7bCi +glC2 + gl Cs + 0{gt) . (421) 

The amplitudes A^, and Cfc depend on the dimension d, and have UV poles at e = 0. 
They also depend on the number of replica n. They are given by the sum of the diagrams 
calculated in section |V} more precisely they are given by 

order k diagrams order k diagrams order k diagrams 

We only need the poles and the finite part of these amplitudes. We find 

Ao =1 , Bo = -n (423) 

n(n—l)l n(n — 1) . 
Ai = 2 e " 2 + ^^'^ ^^^^^ 

n(n-l)(n + 2)l n(n-l)(n-2) , , , 

Bi = -^ ^ - + ^ ^ - + 0{e) (425) 

n{n - l){n - 2){n + 5) 1 n{n - l){-n^ - 6ra + 14) 1 

= 8 ^ + 4 e 

njn - l)(3n^ + (111 - An')n - 210 + Svr^) 
+ ^ + C'(£) (426) 

n{n - l){n + 2){-n^ - 5n + 12) 1 ^ n{n - l){n^ + Qn'^ - 20n + 16) 1 



8 4 e 

n{n - 1) [-3n^ + (-99 + ATT^)n^ - 54n + (456 - 267r2)] 
"I ^ I" ^(^) (427) 

^ 'n,(n — 1) , , 

Ci = ^ ^ ^ (428) 

(-1 + n) n (-7 + 3n + n^) l (-1 + n) n [11 - 10 n + 2r2^ + ln(16)] , ^ ^^^^^ 
(L2 ^ HUl^ej [4:Zy) 

C3 = ^ ^ ' X j (-7 + 3n + n^){-12 + 5n + n')- 

- [-104 + 142n - 58n^ + 3n^ + 2n^ + ln(16) - 8 ln(64) + 2n ln(1024)] ^ 
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+ - 2238 - 2091n + 231^^ + 102^^ + 18n* - MOvr^ + 137nn'^ - 26^^^ - An^n"^ 
6 - 

+ 504 ln(2) - 276n ln(2) - 360 ln(2)2 + 180n \n{2f + 12n^ ln(2)2 + ln(4096) 
+ Oie 



(430) 



Finally, the partition function with one $ operator defined by (|207p is given in the bare 
theory by 



?(2) 
■'3> B 



hgE 



:9i + 0{gt) , 



with 



= 1 

--{n - l){n + 2)^ + (n - 1)(2 - n) + 0{e) 

---{n - l)(n + 2) {n^ + 3n - 7) \ ~ -{n - 1) {2n^ + 5n^ - 24n + 20) - 



e 



+ —{n - 1) (iSn^ - 477^^2 + 96^2 - 777^^ - Qn + 26tt^ - 348) + 0{e) 



(431) 

(432) 
(433) 

(434) 



B. Counterterms and RG functions in the MS scheme 



1. Counterterms 



We use the minimal subtraction (MS) scheme described in sections IV D and IV E The 
one- and two-strand partition functions are renormalized according to (181) and (183) into 
^rHqr! 9r) ^-iid Z^'^\qji ■ ■ ■ ; g^), with the coupling constant g and the momenta q renor- 



malized by ( |173 ) and (180) and the renormalisation factors for the field Z, the coupling 
constant and the boundary Zi of the form (170), (171) and (172). This means that the 



(435) 



renormalized coupling constant is related to the bare coupling constant by 

9b = 9R^g{9R) H9rY^' , 
II being the renormalization mass scale, and that the three amplitudes 

A(^7B)e-2"^^(^-) , B(^e)Z((7H,)e-2"^^(^-) , C((7b) Z(^7h)"^'+'^ e-^"^^(^«) (436) 

considered as series in the renormalized coupling constant (^r, must be UV finite at £ = 0. 
In the minimal subtraction scheme (MS) the counterterms Z, Z^ and Zi are chosen to have 
only poles in 1/e. Putting together ( [419| , ( |42l| ), (|423|-(|430|), ( |435| ) and ( [436] ), we obtain for 
the counterterms 



H9b 

^g{9B 

Zi(5'h 



1+9b 



(n 



+ 9t. 



(n — 1) f A — n 



+ 



n 



1 + 9^^^^ + 9l{n - 2){An - 7) 



3 1 



2e 



+ 0{gl) 
- 0{gl] 



in -I] 



n — 1) f 2n — h n — 2 



9 b: 



+ 0{g[ 



(437) 
(438) 
(439) 
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Similarly, the amplitude 

must be UV finite. This fixes the Z$ counterterms to be 



1 + 



n] 



3(n-2)(n-l) {n-l) 



2e^ 



As 



9l + O [gi) 



(440) 
(441) 



2. RG functions in the MS scheme 



We now use the definition (188) of the (3g function to obtain 

(3g{g^) = -eg^ + [(5 - 2n) + e{n - 1)] gl + (3 - 2n) [(5 - 2n) + e{n - l)]gl + 0{gt,: 



Using (191) and (194) we obtain the anomalous dimension for the field 

7(^7r) = (1 - n)g,, + (1 - n)(3 - 2n)gl + Oigl) , 
and the anomalous dimension for the boundary operator 1;, 



7l(^B 



4 



n) 1 
— S-R H 



Finally, using (203) we obtain the anomalous dimension for the $ operator 

{n - l){2n 



7* (fi-n) = Pgign) ^ ^^^^ = g^{n - I) 



9l- 



11 



dg^ ---^ -' 2 
We remark that we still have at two loop 

^{g^)=A^^{g^)^0{gl) , 
but we do not know if this is a general relation valid at all orders. 



+ 0{gi) 



(442) 
(443) 

(444) 

(445) 
(446) 



C. Counterterms and RG functions in the MS scheme 

We now compute the RG functions in the MS scheme. The counterterms are Z , Zg and 
Zi with g^ = gT^1jgfi~^ and they are chosen to make the three amplitudes in (436) finite at 
e = 0. We obtain 



Z 
Z, 
Zi 



{n-l) 



9b 



(5 - 2n) 



-9r 



{n-A)in-l) (3n-4)(n-l) \ 

2^2 2s fy^^^yy^ 

(2n-5)2 (2n-5)(2n-3) 



+ 



2e 



9l + 0{g-_ 



(n 



11 



(n-l)(2n-5) (n-l)(3n-4) 



9i 



0{gl) 



4e " \ 

We find for the beta function and the anomalous dimension in the MS scheme 
^,(^7r) = -eg^ + (5 - 2n)gl + (5 - 2n)(3 - 2n)gi + 0{gi) 
7(^7r,) = (1 - n)g,, + (1 _ n)(4 - ^n)gl + 0{gl) 

1 



^^7r(1 - ^) + ^^7r(1 - t2)(4 - 3n) + 0((7 



= 9nin - 1) - -(7^(n - l)(4n - 3) + Oigi) 



(447) 
(448) 
(449) 

(450) 
(451) 

(452) 
(453) 
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D. Counter-terms and RG functions in the grand- canonical scheme 



In section VB we have seen that the perturbation expansion in the grand-canonical 
scheme is in terms of the dimensionless combination 



9b ■■= gsiT + nql/Zy 
Using this, we find in the grand-canonical scheme 



+ 9I 



in — l)n 1 , , ^ , , 
{n — 2){n — l)n{n + 5) {n — l)n (n^ + Qn — 14) 

+ - 2)(n - l)n (tt^u + 18n - Stt^ + 102) + 0{e) 



(454) 



0{gi) (455) 



-n + gj, 



[n — l)n 



+ 0{e) 



{n - l)n {n^ + 2n - 6) {n - l)n {2n^ + n - A) 



2^2 



2e 



24 



(n 



- l)n (27rV + 12n2 - lln'^n - 2An + Un^ + 24) + 0{e) 



+ 0{g-_ 



(456) 



C^(fl'B) = 7;9B{n- l)n + g- 



in — l)n {n^ -|- 3n — 7) 
2e 



{n - 1) V + 0{e 



(n - l)n (n^ + 3n - 7) (n^ + 5?2 - 12) 
4^2 

(n - l)n (4n^ + 19^3 - 63^2 + Ihn + 2i 
4i 



+ ^(n - l)n{Ti'^n^ + 90n^ + 2'K'^n^ + 24n^ - 497rV - 102^^ + MOvr^n 



-504n - lievr^ + 192^ + 0{e) 



+ 0{g, 



(457) 



Since this scheme will mostly serve as a check for the integrals in the fixed- length scheme, 
we restrict ourselves to one possible subtraction scheme, namely the grand-canonical MS 
scheme. This gives the RG factors, the index r indicating the grand-canonical scheme. 



Ti 1 

Z- = l + g^^— + gl 



e 

5-2n 

e 



{n-A){n-l) (n-l)(3n-4) 



2^2 



2e 



+ 0{gi) 



+ 9I 



n — 1 2 
^1= 9n^^ + 9^ 



(2n-5)2 (2n-5)(2n-3) 
^2 + Ye 

(n-l)(2n-5) (n-l)(3ra-4) 



8^2 



+ 0{gi) 
+ 0{g 



(458) 
(459) 
(460) 
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The RG functions read 



/9^(^7r) = -9^.e - 9l{2n - 5) + gl{2n - 5){2n - 3) + 0{gt,) (461) 
11, .1 



rig.) = 2 + 2^-(^ - ^) + 2^'^'' ~ ^^^^"^ ~ ^) + ^^^^^ ^^^2) 
7[(^7h.) = \9nil ~n) + ^glil - n)(4 - 3n) + 0(^?^) . (463) 



They are identical to (450)-( [453 ), confirming our calculations in the fixed-length scheme. 



E. Fixed point and critical exponents 

We now analyze the RG fiow within an e expansion as in LW [2Hj and in our letter 
First, we give all results at arbitrary n, before specifying to = in the next subsection. 
For small e > and n > —5/2, thus especially for the limit of interest n —>■ 0, there is an 
IR unstable fixed point g*. In the MS scheme this fixed point is at 

This freezing transition separates the molten (weak coupling) phase with < g < g*, where 
disorder is irrelevant, from the (strong coupling) glass phase with g > g*, where disorder is 
strongly relevant. The anomalous dimensions at the fixed point give the critical exponents 
at the glass (freezing) transition. 



The correlation-length exponent u defined in eq. ( 14 ) is 



1 1 2n — 3 , , , , 

The scaling dimension of the overlap operator, which gives the decay of the pair correlator, 
is 

1 2^ — 3 

e*=2-^'^{g*) = 2-- = 2-e + ^-^ + 0{e') . (466) 

It is in agreement with our result announced in [H] (note that contains a sign typo in 
the intermediate identity of Eq. (20), but the final result at order e"^ for 6* is correct). The 
anomalous dimensions of the and contact operator $ are respectively 

This gives the scaling dimension for r, 

^ -2~y-2 + 4;7^"+2(5-2n)2" + > ^^^^> 
and the dimension of the contact operator $, 

^ ^ + 2^^* ^10-4n 2(5-2n)2^ ^ ■ ^ ' 
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Finally, our explicit calculation shows that the scaling relation (212), (* = {2 — b)x*, holds 
at two loops. Indeed, 

^* + p* _ (2 _ e)x* +p* = 2 + 0{e') , (471) 
in agreement with the exact scaling relation ([8|. 



F. Results for the disordered system {n = 0) 



The physically relevant case is the zero-replica limit n = 0. There we find 

* 1 3 
u* = 

e 5 



^-1 + 0{e) (472) 



9* = 2-e--e^ + Oie^) (473) 
5 

P* = 1 + + + 0{e') (474) 

X*^^ + ^e + l.^ + 0(.3) (475) 

C = 1 - - + 0{e') . (476) 



The physical case for the RNA folding model corresponds to e = 1. 



G. Physical interpretation 



In fig. [19} we have plotted, for n = 2, n = 1 and n = 0, 9* and p*, as obtained through 
our 2-loop calculation (solid lines), as well as the resummed 1-loop results of [28] (dashed 
lines) 



1 + 1 + ^^^ 



~'~ 5-2n ~'~ 5-2n 

The latter are interesting, since they are exact for n = 2, where they can be obtained by an 
exact solution of the corresponding self-consistency equations [2H1 SH]- For n ^ 2, they are 
expected to be better than the pure 1-loop results, which are obtained by extrapolating 9* 
and p* linearly around e = 0. 



Consider now fig. 19 for n = 2. The dashed lines (resummed 1-loop results) are exact. 
However since p* > po = 1 + | (dotted line), the exponent p* gives only the subleading 
corrections to the contact-probability, and the dominant contribution comes from the mean 



field exponent po = 1 + e/2 (represented by the dotted line on fig. 19, first graph), instead of 
p*. Since base-pairings common to two replicas are a subset of base-pairings in one replica, 
the pairing probabilities satisfy 

(^«/3(s,t)) < (0«(s,t)) (478) 

and as a consequence 

9* >p* . (479) 



Note that 9* > po, consistent with (479), and the inequality is attained for 6=1. One can 
show [in] that for £ > 1 both exponents take the value |, and that there are logarithmic 
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n = 



FIG. 19: Results for 9* (black) and p* (grey) at 1- (dashed) and 2-loop (solid) order, for n = 2, 
n = 1, and n = 0. The dotted line for n = 2 is the value of po, which is an upper bound for p* . 
For interpretation see main text. 



corrections at e = 1. We also see that the 2-loop corrections are rather large, and a resum- 
mation appropriate. However, as we will see below, for the smaller n we are interested in, 
this will not be necessary. 

Consider now n = 1. This case is peculiar, since the disorder lives on pairs of distinct 
replicas, thus "does not exist" for = 1, as can be seen e.g. from the factor of {n — 1) in the 
contribution to the anomalous dimension 7$ in eq. (468). As a consequence, p* = po exactly, 
which is respected by our perturbation theory. The (3 function can nevertheless be defined. 



as the limit of n — * 1, resulting in the exponent 9* plotted in fig. 19 In order to respect 



(479), one has to set 9* to p* = po as soon as the two curves cross, i.e. for e > Ec = 0.589. 
We call Ec the upper critical dimension, i.e. the dimension at which at the freezing-transition 
the two replicas get locked 

We now arrive at the physically relevant case n = 0. The first observation is that p* < po, 
thus disorder makes the pairing probability more long-range correlated, and p* is the leading 
exponent for the latter (this is in fact true as soon as n > 1). As in the case n = 1, beyond 
the upper critical dimension, i.e. e > Ec = 0.592, the inequality (479) is violated. Following 
[28] , we assume that the value of p* can again be extrapolated analytically beyond Ec, and 
that 9* is given by p*. We have no proof that this assumption is correct, but expect that 
due to the very character of the e-expansion, which states that the topology of fixed points 
and their corresponding exponents evolves smoothly upon a change in parameters (here n 
from n = 1 to n = 0), the possible error is at least small. Last not least, we remark that 
2-loop corrections (especially for p*) are small, thus no resummation is necessary. These 
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arguments lead to the following value for 6* and p* at the transition 

e* = p* = 1.36 . (480) 

Our discussion above implies that at the transition the two replicas are already locked 
together, i.e. for a given disorder one possible fold dominates the partition function. This is 
not what we would expect. Indeed, at high temperatures, in the so-called molten phase, all 
possible folds are equally probable; and one expects that lowering the temperature, at the 
transition only a subset remains relevant, whereas at zero temperature a single configuration 
dominates. This naturally leads to three different sets of critical exponents: above, at, and 
below the transition. However, since our theory suggests that at the transition as in the low- 
temperature (glass) phase the same single configuration dominates the partition function, 
we are led to conjecture that in the glass phase 

pglass _ pglass _ ^* _ ^* _ ^4g^^ 

This "locking scenario" was first proposed by Lassig and Wiese p5|. Our final result for p 
and ( at the transition and in the glass phase is 

Pglass = P* = 1.37 ± 0.01, Cgiass = 2- p* = 0.63 T 0.01 , (482) 

where the central value and the reported error is an estimate based on the three possible 
Fades for C*{^)y the Fades (2,0), (1,1), and (0,2). The error due to the neglect of higher- 
order corrections is difficult to estimate, and not given here. Numerical results obtained by 
Krzakala et al. [H] in agreement with Bundschuh et al. pT| [55] give 

Pglass ^ 1.34 ± 0.003, Cgiass ^ 0.67 ± 0.02 . (483) 

These numerical results compare favorably well with those from the renormalization group, 
using the locking argument. 



VII. RANDOM RNA UNDER TENSION 

An interesting question is what happens when a RNA molecule is pulled at its both 
ends. This problem has been studied in flSl [261 [27] both numerically and analytically in 
the molten phase, and experimentally in [22]. In |15] a generalization of the field theory was 
proposed for random RNA under tension, and the force-induced denaturation transition. In 
this model the pulling force / is treated as a small perturbation. A RG calculation to first 
order in the force was performed to study the denaturation transition and its interplay with 
the disorder-induced freezing transition. 

In this section we recall the definition of the model, show that it is renormalizable for 
e = as for the tension-less random RNA model, and give the details of the calculation of 
the RG flow and of the critical exponents at two loops. For a more detailed justification we 
refer to |15], where some of the results were already announced. 



The model 



The discrete model for random RNA under tension is obtained by adding to the disorder 



Hamiltonian 7i[$] given by (50) a new term proportional to the external force / of the form 

Hf_[<l>] = -/5^5^A,(z) . (484) 
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► r 



FIG. 20: (a) An open planar RNA structure, (b) arch diagram and (c) corresponding height rehef. 
The free part of the strand is in red, corresponding to height /i = 0. 



Here, the "free-strand operator" Aa{i) is defined for each replica a and for each site i as 

A„(2) = <^ ^ ' (485) 

I otherwise , 

where the height h{i) of the site i is defined via the height representation (|6]) of the planar 
pairing system, 

j<i k>i 

The graphical interpretation of this operator is depicted on fig. 20 

In the continuum limit, this operator corresponds in the planar random- walk representa- 
tion to the 1-point contact operator 

a,b a,b 

which represents the interaction of the random walk r(s) with an impurity fixed at the 
origin r = 0. The dressing by the auxiliary fields 7, 7 (which are the same fields as in 
the previous sections) is introduced in order to eliminate non-planar diagrams. The force- 



insertion operator A.a{u) is represented diagrammatically on fig. 21 by a single dashed line 
connecting the random walk (represented by the full line) to the impurity (represented by a 
cross). The planar dressing by the auxiliary fields is represented as in the force-free case by 
additional lines carrying color indices a (a = 1, ... , A^). Their role is to replace the lines by 
fat lines so that the planarity constraint is implemented in the limit of ^ cxd. 

This planarity constraint is required since the force term only acts on the free part of the 
RNA strand. It is easy to see that once we include the disorder induced interaction between 



replica, non-planar diagrams such as those of type (b) in fig. 22 represent force terms acting 



on the non-free part of the RNA strand, since they are inside an interaction arch. Only 



planar diagrams such as those of type (a) in fig. 22 are to be taken into account. 



The force is thus taken into account in the action by the additional interaction term 



A^{u) (488) 

0<M<L 



81 



(7 X 



FIG. 21: Diagrammatic representation of the dressed force insertion operator Aa{u) 



X 



X 



(a) 



(b) 



FIG. 22: (a) A force insertion on the free strand is represented by a planar diagram, of order 
0{1). (b) A non-planar diagram represents a force insertion on a non-free part of the strand and 
is of order 0{N-^). 



In the RW picture it corresponds to the addition of an attractive short-ranged pinning po- 
tential at the origin r = for each polymer. The potential is attractive because of the minus 
sign in the action (the pulling force is /o > 0). 

The full action of the model with disorder and force is 



S — So -\- (Sint + 



(489) 



and depends now on two coupling constants, the bare disorder-induced coupling constant 
go > 0, and the bare pulling force /o > 0. The naive (engineering) dimension (in units of L) 
of the force-insertion operator Aq,(w) and of the bare force /o are respectively 



A/- 2 



(490) 



The force / is thus relevant if £ > 0, exactly as the disorder couphng go. 

Finally we consider observables in the force model. The length of the free open strand in 
the discrete model (for a single replica) is 



(491) 



Therefore in the continuum model it is 



'-'free 



= / A(«) , 

Jo<u<L 



(492) 



and it has engineering dimension e/2. Secondly we consider partition functions for open 
strands, that we describe in the next section. 
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B. Perturbation theory 

We can now build a systematic perturbation theory by expanding observables as a power 
series in both qq (the disorder strength) and /o (the pulhng force). The diagramatics involves 



both the two-replica bi-local vertex of fig. 10, representing the overlap operator "^a/s in 



■^forcc • 



and the single-replica force-insertion vertex of fig. 21, representing the operator in St, 
This last term breaks the translational symmetry of the / = model, since it corresponds 
to the addition of an attractive short-ranged pinning potential at the origin r = for each 
polymer. Therefore one must be careful when dealing with the zero modes in perturbation 
theory. 



1. One strand partition function: Zero-mode decomposition 

We consider the partition function for a single (bundle of n replicas of an) open strand 
with length L and external momenta q", qf for each replica a. It is defined as 

S(i)(qi,q2;qo,/o) = / V^] e"^!''-! e 'e^^' (493) 



In the perturbative expansion of (493) in powers of go and /o, we must be careful how 
many different replicas get at least one force insertion, since as soon as a replica gets a 
force insertion vertex, no 5'^(qi + q2) occurs (since there are no translational zero modes 
any more). Therefore we decompose the partition function S(i)(qi, q2; qo, /o) in a sum over 
sectors labeled by the number k of different replicas with at least one force insertion, and 
factor out for the n — k remaining replicas the 5(qi+q2) factor coming out of the translational 
zero mode, so that < k < n: 

n 

2(i)(qi,q2;^o,/o) = Yl ^'^(qi + ^2)^""''^ S(i,fc)(qi,q2;^o,/o) • (494) 

k=0 

Each sector (k) gives a finite partition function S(i fc)(qi, q2; (^o? /o) and the translational 
zero modes are taken into account in the factors of 5'^(qi + q2). 

One should keep in mind that, as in the force-free case discussed in the previous sections, 
each replica may carry different incoming external momenta (q", q2 ); for simplicity we take 
these momenta independent of the replica index a. This is why we must keep the dependence 
on the two momenta qi, q2 for '^(i^k)- Note also that for the sector k = (no force insertions) 
we recover the partition function for the single free strand (no force) 

H(i,fc=o)(qi, -qi; ^o, /o) = ^(i)((qi; ^o) • (495) 

As we shall see, only the sector {k = 1) is in fact needed to study the renormalisation of the 
model. It is of course possible to define and study multi-strand partition functions in the 
presence of a force, and to treat properly the zero modes. Fortunately this is not necessary 
either. 



2. Diagrammatics and UV divergences 

We now compute the partition function S(i fc) up to second order in the couplings go 
(disorder) and /o (force). The corresponding diagrams are detailed below. Their amplitudes 
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are calculated in dimensional regularization using the same methods as in the force-free case. 
As discussed above, the force-insertion vertex is associated with a delta function in position 
space, which is represented by an external momentum-vertex insertion 



d/2 



^ gikarQ,(?i) 



(496) 



As explained in [45J, using this representation for the force insertion vertex, together with 
that for the bilocal contact operator 



TC 



d/2 



(497) 



taking the e.v. w.r.t. the free field r, and then integrating over the momenta k, we obtain an 
integral representation for the amplitudes of the diagrams in terms of the relative positions 
of vertices (the m's and f 's) on the line [0, L]. 

Let us look at a few simple cases to see how this works and where the UV divergences 



occur. The complete list of diagrams is given in section VII D 2 

The diagrams without disorder {g = 0) and {k = 1) (force insertions on one replica) 
are the diagrams for the theory of a free polymer pinned by a single impurity. They are 
represented with the corresponding amplitude, qi and q2 are the momenta entering at the 



both ends, and we drop the combinatorial factors (given in section VII D 2 ) . The points are 
labeled x, y and z (if appropriate) from left to right. 







^ > 

i 
1 




> 


c > 


c > 





e ^1 e 



0<x<L 



e ^1 e 



\y - x\ 



0<x<y<L 



' 0<x<y<z<L 

The diagrams of order f x g are with the same convention for the points x, y and z: 



(498) 
(499) 
(500) 



-<if^p-<i2iL-^+y-^)p-Qhp-'i2iL-^)\ 



e ■'^"■e ~ "^'e ■'^"e ~'\z-y\ ^ (gi + ^2) 



0<x<y<z<L 
>(. 



(501) 



^ ^ 



0<x<L 



0<y<z<L 



(502) 
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(Here y and z are the points on the first two rephcas, and x is the point on the third rephca.) 
These amphtudes can be computed exphcitly, and contain short-distance UV divergences 



(UV pole) at e = 0. For instance (499) contains a pole at £ = coming from j/ — > which 



is proportional to the amplitude of (498). Therefore it is proportional to the insertion of a 



force operator A(x), and can be absorbed into a renormalisation of the force /. 



Similarly, (501) contains (in addition to the poles coming from z ^ y and {y, z} 



L), a pole coming from {y, z} x, which is proportional to (498), therefore also to a 
renormalisation of /. 

These new UV divergences can be analyzed by the Multilocal Operator Product Ex- 
pansion techniques already used for the tension-free model. One must consider the short- 
distance behavior of products involving both replica-interaction operators \E'a/3(M, f) and 
force-insertion operators Ao(u). We do not discuss this MOPE in more detail, but mention 
that at first order, in addition to the standard MOPE, 

^xr^^l+r^H , ^x^^l+r^H , ^I'x^^^H , ^'xlfo^l^H (503) 

the following additional terms appear: 

*xA->A + -- -, AxA^A. (504) 

This implies that there is a force renormalisation induced both by the force / and by the 
disorder g, but that the force / does not renormalise the disorder g. 

C. Renormalization 

1. Renormalized action and beta functions 

The renormalized action for the model of open RNA strands under tension is thus (we 
omit the auxiliary fields 7 and 7) 

^n[^-] = Y. [h^-Y.9^^9f^'' ff $^, + $^2Zi-5^/,Z,/.-/2 /■ A^ (505) 

a ^ a<f3 •' ^ a a ^ 

The counterterms Z = Z(5fR), Z^ = ZglgT^) and Zi = Zi(5fR) are the same as those for the 
force-free model, since they are not changed when / > 0. /a is the renormalized force and 
Zj = Z/(/r, (7r) is the new force-renormalisation factor. It depends both on and on g^^. 
For the disorder-free model g^i = 0, one recovers the renormalized theory for the model of a 
random walk pinned by an impurity [291 EHl IMl EH [56j . 

The renormalized action can be written as a bare action iSnfra] = iSB[rB], with the same 
bare fields and bare couplings rg and g^^ as before, but with a bare force /b given by 

g^ = (7rZ,ZV-' , /b = /rZ/Z'^/ (506) 

The renormalized functions fc)(qi, gn, /r) calculated with the renormalized action are 
UV finite. They are given in terms of the bare functions fc)(qi5 q2; ds, /b) (calculated with 
the bare action) by 

Sfi,.)(qi,q2;^7K,/K) = Z-'^/'e-'-^^ Ef,^,^iZ-'/'ci„Z-'/\,; g,, (507) 
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The factor comes from the zero modes. 

The RG beta function for the disorder couphng (yfj^, f3g{g^), is unchanged. In addition 
there is a RG beta function for the force, which is defined as 



dfn 



dfj, 



9b Jb 



/r 
/r 



e d d\nZ dlnZf 



{2 + e)Xr{9n) - Pgign) 



rflnZ 



(508) 



2. MS versus MS schemes 



As for the force-free model, we consider the two different subtraction schemes MS and 
MS. In the MS scheme, the counterterms Z, Z^, Zi and Z/ are chosen such that they contain 
only pure poles in e, and no finite or analytic term at £ = 0. It is more convenient to consider 
the MS scheme, where the coupling-constant renormalisation factor Z^ = "Lgl?^^ has only 
pure poles. This scheme is easily extended to RNA under tension: 



Z/ = Z /Z^+"/2 g^^^j^ ^^^^ f^^f^Zf fj,-'/^ 



(509) 



In the MS scheme, it is Zj, together with Z^, which has only pure poles and no finite analytic 
part at £ = 0. The /3j function becomes 



5/r 



d/j, 



fn 



9b Jb 



e ^ , , dlnZf 



(510) 



D. Two-loop calculations 



1. Principle 



Since we have already computed the boundary and the wavefunction counterterms (they 
are independent of /o), it is enough to compute the bare ^^(qi, q2; (^b, /b) function for 
k = 1 (force insertions on a single replica) and at zero momenta qi = q2 = 0, and to isolate 
its UV poles at £ = 0, in order to compute the force counterterm Z/ and the RG functions 
for the force. 

Further simplifications occur, since the contribution of the diagrams with a single force in- 
sertion are related to the contributions of the corresponding diagrams with no force insertion 
at first orders in q, which we have already calculated. For instance 



+ 



ld_ 
'2d^^ 



(511) 



5=0 



q=0 



Similar identities hold for the two-loop diagrams. 
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2. Diagrams 

With these remarks we now collate the diagrams Qi, their symmetry factors c{Qi) and 
their amplitudes Gj for the bare function Si fe=i(qi — q2 = 0; 5(3, /b) (with the strand length 
set to L = 1) up to two loops 

5(i,fc=i)(qi = q2 = 0;^,/)= J] 9^f^c{g,)G,. (512) 

diagrams Qi 

When possible, we express the new amplitudes Gi in terms of the already calculated am- 
plitudes Ai and Bi of the force-free model. We do not give results for the grand-canonical 
scheme. 

Order /q: 



c(^i) = n , Gi = 1 



(513) 



Order /q: 



52^ 



Order fogo- 



C(^2) = , 6*2 = / (1 - U)U 

J u 



+ 



17 



c(^?3) = n(n - 1) , Ga = ^ - -S2 



-(i/2 



{2-d/2){l-d/2) 



n{n — l){n — 2) 



(514) 



(515) 

Gi = ^1^2 (516) 



Order fl 



XXX 



c(^5) = n , G5 



r(i-d/2)- 
r(4 - d) 



(517) 



Order x qq. 



/66 



/6c 



, ciSea) = n{n - 1) , G^a =^ / r-''s-''/'{l 

^ Jr+s+t=l 

r(-i-£)r(-5/2) 



r — s 



r(2 - f ) 



6a 



:7 



:7 



, c(e?6c) = n(n - 1) , Gee = / 

-'t+M) + X=l 

r(i-d)r(2-rf/2) 



r(5 - 3d/2) 



(518) 
(519) 



(520) 



Order fo x g^. 



n{n — l)(n — 2) 



G7 = 6-2^2 (521) 



+ 



+ 



ciGs) = n(n - 1) , Gg = A3 - S3/2 

e?9 = 



+ 



:7 



c(6?9) =n(n-l) , 6-9 = ^4-54/2 

X 



+ 



+ 



+ 



+ 



+ 



7 



(522) 



(523) 



, ciQio) = n{n - 1) (n - 2) , do = 3A, - B, (524) 
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+ 



+ 



(525) 



+ 



+ 



+ 



^" ,7'-, ) 



ciQii) = n{n - 1) (n - 2) , Gn = SAg - 

X 

Gl2 = \ 



(526) 



^13 = 



Ql4 = 



'15 



^16 = 



Gn = 



^/^ X n(n- l)(n-2)(n-3) _ j 

, C(C/i2j = , Gi2 = A2G3 



(527) 



, c(gi3) = "^"-^j^"-^\ G^s-Gaa (528) 



, c(gi4) = , Gu^G.A, (529) 



, c(^?i5) =n(n-l)(n-2)(n-3) , = GiAg 

(530) 

, c{g,e) = n{n - l){n - 2){n - 3) , = Gi^ 

(531) 

n(n- l)(n-2)(n-3)(n-4) 



, c(^i7) = 



:7 



Gir =GiAy 
Gathering the results we get 

S(i,fe=i)(qi=q2=0; g^, /b) = Eo,i/b + Ei,i^b/b + ^o,2fi + Ei,2^b/b + E2,i^^/b + ^0,3/! 



(532) 



(533) 
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with (keeping the singular and finite part in e for the amphtudes E. 

2 



Eo,i = l , 

+ n - 2 1 



Eo,2 = -- - 1 



4 4 ^ 

Eo,3 = ^ + -- ^+ 4 + 



Ei,i 

E2,l 



{n - 2)(n - 1) + 



2e 2 

(n - l)(n + 2)(n(n + 5) - 12) (n - l){n{n{n + 6) - 20) + 16) 

+ - 1) (3n(n(3n + 31) - 10) + 7r2(26 - n(4n + 7)) - 312) + • • 

in + 2)in-l) in-A)in-l) I, , 
Ei,2 = + ^ - 7 3n - 27r2 + 33 (n - 1) + 



(534) 
(535) 

(536) 
(537) 



3. Beta functions 



Starting from the relation (507) between the bare and renormalized functions S(i ^=1), 
using the already obtained two- loop expression for the counterterms Z, Zi and Zg, we are 
left with the determination of the force-renormalisation factor Zj. This counterterm is 
determined by enforcing that the renormalized function j^^-^-^ has no poles in e. In the 

MS scheme this yields 



Z/(^a,/a) = l + /R- + /a^ + /R^7, 



1 — n 



+ 9i 



[n- l)(2n-5) 

45 



(538) 



We find for the 2-loop beta function for the force /r in the MS subtraction scheme, 



/3/(^7k,/k) = -;^/r + /r,+ 



2 I 3(1 -n) (i-^)(5- 2n)^ , 



(539) 



We recall that the beta function for the disorder coupling g^^ in this scheme is 

= -eg^ + (5 - 2n)gl + (5 - 2n)(3 - 2n)gl + ■ ■ ■ . (540) 



E. Freezing and denaturation transitions 

We now discuss the physical applications, following [35]. To apply our calculation to the 
problem of random RNA under tension, we must take the n = limit, and set e = 1. The 
RG functions become at n = 



PfignJa) 



-eg^, + 5gl + 15gl + 



-Ml + 



(541) 
(542) 



Note that there is no term in the MS scheme. 



We recall, that equations (541) and (542) give the RG fiow, when going to large scales 



(large RNA molecules). For small g^^ and /r, the fiow is dominated by fixed points of order 
S'r ~ /r ~ ^) and can thus be studied in a Wilson-Fischer small e-expansion. This fiow 
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FIG. 23: RG-flow as given by equations (542) and (541), with fixed points (red) and separatrices 
(thick orange/grey hnes). 



in the ((/r, /r) plane is depicted on figure 23 Besides the unstable Gaussian fixed point 
O = (0,0), there are three non-trivial fixed points, which read at order e^: 

e ?>e^ 

freezing transition f.p. F fi'p = , /f = (543) 

denaturation transition f.p. D = , /d = ^ (544) 

bicritical f.p. B ^*=--— , /* = ___ (545) 

D is the fixed point for the tension-induced denaturation transition of homogeneous RNA 
(i.e. for homopolymers) . F is the fixed point for the force-free RNA freezing-transition. B 
is the new unstable fixed point of [13], corresponding to a bicritical freezing+denaturation 
point. 

The denaturation exponent 7 at weak disorder is given by the derivative w.r.t. /r of 
Pf{gn, fa) at the fixed point D, 



7 = 



= = 0.5 , (546) 

/ e=l 

D 



(this is in fact an exact result). At the bicritical point B this denaturation exponent is 
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modified to 



7 



dpfigK,fR 



+ 



10 



0.6 



(547) 



Assuming the locking scenario of Lassig-Wiese, this exponent should be equal to the critical 
exponent for the tension-induced denaturation transition of RNA in the frozen phase {g 
large). Our 2-loop estimate 72_ioop — 0-6 is to be compared with the result of numerical 
simulations of Krzakala et al. [13], 7sim — 0.7, and is definitely better than the 1-loop 



estimate 



loop 



0.5. 



VIII. CONCLUSIONS AND PERSPECTIVES 

In this article, we have considered the folding of RNA strands with random pairing ener- 
gies, as a model for folding of random RNA sequences. We have established the existence of 
a phase transition from a high-temperature/low-disorder phase with a pair-contact exponent 
p = I to an exponent of p* ~ 1.36 at the transition. This was achieved via a designed field 
theory, modeled along the lines of self-avoiding polymers and membranes. It allowed us to 
use the tools of the multi-local operator product expansion, show renormalizability, and to 
obtain results at 2-loop order. 

Our considerations here are for RNA-molecules in thermal equilibrium. An important 
question is whether the dynamics of RNA folding [21 157] , especially after a quench can be 
modeled too, or melting scenarios a la [58j explored. 

The order of the phase transition may depend delicately on the role of excluded-volume 
interactions. It is now possible to ask questions about the tertiary structure, using as a 
starting point our field theory and the ensuing statistics of branching, or equivalently the 
statistics of the RNA fold seen as a tree. 

Our results suggests that steric interactions should be relevant. Indeed, the internal 
fractal dimension is cif = l/( ~ 1.56. If the embedding of the tree in external space is 
dominated by entropic effects (mean field), so that the tree forms a "blob", the fractal 
dimension of the tree in bulk space should be duoh = 2(if ^ 3.12. Steric interactions are a 
priori relevant if the dimension of the embeding space is demh < 2(ibiob, thus are important 
for long strands, even though their effect may not be pronounced since 2(ibiob — c^emb is small. 



The numerical results in eq. (483) are indeed consistent with 2(ibiob — c^emb = 0, which would 
lead to logarithmic corrections only. 

Finally, in view of the results of [2] it will be necessary to reanalyze the role of knots and 
pseudo-knots, even though their existence may not change the asymptotic scaling. Our field 
theory is capable of achieving this, by including corrections in 1/A^, i.e. corrections to the 
planar limit, systematically, in the spirit of what was done in [6] for homopolymers. 
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